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Abstract 

Neutralino dark matter, and in particular different aspects of its detection at neutrino 
telescopes, has been studied within the Minimal Supersymmetric extension of the Standard 
Model, the MSSM. 

The relic density of neutralinos has been calculated using sophisticated routines for 
integrating the annihilation cross section and the Boltzmann equation. As a new element, 
so called coannihilation processes between the lightest neutralino and the heavier neutralinos 
and charginos have also been included for any neutralino mass and composition. 

The detection rates at neutrino telescopes have been evaluated for neutralino annihilation 
in both the Sun and the Earth using detailed Monte Carlo simulations of the whole chain 
of processes from the neutralino annihilation products in the core of the Sun or the Earth 
to detectable muons at a neutrino telescope. 

A comparison with other searches for supersymmetry at accelerators and direct dark 
matter searches is also given. 

The signal muon fluxes that current and future neutrino telescopes can probe and the 
improvement in sensitivity that can be achieved with angular and/or energy resolution of 
the neutrino-induced muons has also been investigated. 

The question of whether the neutralino mass can be extracted from the width of the 
muon angular distribution, if a signal flux is observed, has also been addressed. 
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Chapter 1 

Introduction 



Many cosmological observations show a definite need of dark matter, which can make up 
more than 95% of the mass in the Universe. One usually defines = p/pait where p is 
the density in the Universe and p cr i t is the so called critical density for which the Universe 
would be flat. Rotation curves of galaxies indicate that 

ft > 0.1 (1.1) 

in contrast to the luminous mass density 

ftluminous < 0.01 (1.2) 

which clearly indicates the existence of dark matter. Moreover, motions of galaxies in clusters 
and superclusters indicate that 

ft> 0.2-0.3 (1.3) 

or maybe even higher. ft is also bounded from above due to the age of the Universe being 
at least 10 10 years f|, 

ft/i 2 <1 (1.4) 

where h is the Hubble constant in units of 100 km Mpc -1 s _1 . The value of h is still a bit 
uncertain with different experimental determinations ranging from 0.5 to 0.9. 

What can this dark matter be then? Studies of Big Bang nucleosynthesis predict values 
of the abundances of 2 H, 3 He, 4 He and 7 Li which, when compared with observed abundances, 
give § 

0.008 < fl hmyon h 2 < 0.024. (1.5) 

We thus see that the dark matter cannot be made up of baryons, but some more 'exotic' relic 
from the big bang is needed to explain the high values of ft observed. There exist several 
hypothetical candidates of which a Weakly Interacting Massive Particle (WIMP) is a major 
candidate. WIMPs will freeze out in the early Universe when they are non-relativistic and 
their relic density is approximately given by U 

j2 3 x 10" 27 cm 3 s- 1 . . 

"wiMpn. — -, r (1.6) 

where (<ja v ) 1S the thermally averaged annihilation cross section. A weakly interacting parti- 
cle is expected to have an annihilation cross section of the order of (ctav) ~ a 2 w / (100 GeV) 2 ~ 
10~ 25 cm 3 s _1 which gives an ft of about the magnitude wanted. Hence, if there are any 
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WIMPs left from the big bang, they are expected to have a relic density that can be enough 
to explain the dark matter problem. 

One of the leading WIMP candidates is the lightest supersymmetric particle, the neu- 
tralino. Assuming that the dark matter is really constituted by WIMPs (or more specifically 
by neutralinos), how can one find them? First of all, if WIMPs constitute the dark matter, 
they will have clumped together making up a halo of the galaxy (containing most of our 
galaxy's mass). Hence they will be all around us and the Earth (and the solar system itself) 
will move through this halo during its motion through the galaxy. 

There are in principle two different kinds of experiments proposed to search for the dark 
matter in the halo of the galaxy: direct and indirect searches. In direct experiments one 
looks for these WIMPs passing by a detector and scattering off some nucleus. This scattering 
can be detected and, if found, would be an evidence for WIMPs in the galactic halo. In 
the indirect searches, one looks not for the WIMPs directly, but for signals coming from 
annihilation of two WIMPs. For example, their annihilations in the halo will result in a 7- 
and p-flux which can be searched for. These WIMPs can also get elastically scattered while 
passing the Sun or the Earth and get gravitationally trapped. They will then accumulate 
at the center of the Sun and the Earth where their annihilation eventually will produce 
neutrinos which can be detected. This last indirect way of searching for WIMPs is the main 
topic of this thesis. Though most discussions will be devoted to neutralinos, many of the 
results in this thesis will be applicable to any WIMP. 

In Chapter || the Minimal Supersymmetric extension of the Standard Model (MSSM), 
in which we work, will be defined, in Chapter || the experimental constraints on the MSSM 
will be reviewed, in Chapter || a detailed calculation of the relic density of neutralinos will 
be performed and in Chapter |^ the expected neutrino-induced muon fluxes at neutrino 
telescopes will be evaluated. We then close by some concluding remarks in Chapter 0. 



Chapter 2 

Definition of the MSSM 



2.1 Introduction 

We work in the Minimal Supersymmetric extension of the Standard Model (MSSM) with 
N = 1 supersymmetry generators and we will essentially follow the notation of Ref. 
We will only give a short introduction to supersymmetry phenomenology in this chapter and 
the interested reader is referred to Ref. for more details. 

Supersymmetry is a symmetry relating fermions to bosons such that for each fermionic 
degree of freedom there is a bosonic degree of freedom. This extends the particle content of 
the Standard Model (SM) such that each particle in the SM has a corresponding superpartner 
(or partners). More specifically, the particle content in the MSSM is the same as of the SM 
plus the superpartners and two Higgs doublets (instead of one as in the SM). Two Higgs 
doublets are needed to give mass to both up- and down-type quarks and will result in five 
physical Higgs bosons. If supersymmetry were unbroken, a SM particle and its superpartner 
would have the same mass and quantum numbers (except for spin). Since we haven't seen 
these particles, we can conclude that supersymmetry is broken at the energies probed by 
present accelerators. 



In Table 2.1 we list the 'normal' particles and their corresponding superpartners. Note 
that some 'normal' particles have more than one superpartner, e.g. each quark has two 
squarks, qx and qn as superpartners, but the number of degrees of freedom (2 for the quark 
(spin i) and 1 for each squark (spin 0)) sums up to be the same for the normal particle 
and its superpartner(s). The general notation is to have a tilde on the symbol for the 
superpartners, but for the charginos and neutralinos we will usually drop the tilde since 
there is no risk for misinterpretations anyway. 



2.2 The superpotential, supersymmetry breaking and 
imparity 

To write down the Lagrangian for the MSSM, one should introduce the superficld formalism. 
This is not within the scope of this thesis and we will just write down the superpotential 
and the soft supersymmetry breaking potential for reference and the reader is referred to 
Ref. [§U for details. 

The superpotential is given by 

W = eij (-e R Y E V L H{ - d* R Y D ^ L H{ + u^Y a qi#| - ^H{W 2 ) (2.1) 

7 
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Normal particles/fields 



Supersymmctric partners 
Interaction eigenstatcs 



Mass eigenstates 



Symbol 


Name 


Symbol 


Name 


Symbol 


Name 


q = d,c, b, u, s, t 


quark 


QL, §R 


squark 


qi > 92 


squark 


I = e, fi,r 


lepton 


II, Ir 


slcpton 


h, h 


slepton 


v = v & , , u T 


I1L. U. LI IIILl 




olIL LI Lrl 11IU 




bile LI LI ILIO 


a 


gluon 


9 


gluino 


9 


gluino 


w ± 


W^-boson 


w± 


wino 






H- 


Higgs boson 


#r 


higgsino 


> xf, 2 


chargino 


H+ 


Higgs boson 




higgsino 






B 


B-field 


B 


bino 






W 3 


VK 3 -field 


w 3 


wino 






H° 


Higgs boson 
Higgs boson 
Higgs boson 


m 


higgsino 


> X?,2,3,4 


neutralino 




m 


higgsino 







Table 2.1: The 'normal' particles and their superpartners in the MSSM. 



where i and j are SU(2) indices, the Yukawa couplings Y are matrices in generation space 
and e, 1, u, d and q are the superfields of the leptons and sleptons and of the quarks and 
squarks. The lefthanded components are SU(2) doublets and the righthanded are SU(2) 
singlets. 

We then introduce all possible soft supersymmetry breaking terms (without violating 
gauge-invariance or breaking baryon or lepton number) in the potential 



V soft = ey (e R A E Y E ¥ L H{ + d* R A D Y D q^H( 
-vkAuYu&Hi - B[iH\H 3 2 + h.c.) 



+H[*mlH\ + H%m\n\ 

+qi*M|ql + l»M|li + u^IV^u fl + d^M 2 D d K + e^M|e fl 

+ ^M 1 BB + ^M 2 (W 3 W 3 + 2W + W~^j + ^M 3 gg (2.2) 

where the soft trilinear couplings A and the soft sfermion masses M are matrices in gen- 
eration space and the fields e, 1, u, d and q are the scalar components of the superfields 
corresponding to the superpartners of the SM. The L and R subscripts on the sfermion 
fields refers to the chirality of the fermion they are superpartners of. B, W 3 and W ± are 
the fermionic superpartners of the SU(2) gauge fields and g is the gluino field, yu is the 
higgsino mass parameter and Mi, M 2 and M 3 are the gaugino mass parameters. B is the 
soft bilinear coupling and mi ,2 are Higgs mass parameters. Notice that we have now intro- 
duced many new parameters, but this is the price to pay until we know how supersymmetry 
breaking occurs. 

The superpotential and soft supersymmetry breaking potential we have now introduced 
are not the most general ones, unless we assume that the so called i?-parity is conserved, im- 
parity is a discrete symmetry being 1 for 'normal' particles and —1 for their superpartners. 
i?-parity has to be put in by hand and if conserved automatically prevents baryon and 
lepton number violation which would otherwise be allowed unsuppressed at tree level, im- 
parity conservation also implies that the Lightest Supersymmetric Particle, the LSP, is stable 
which is a very welcome consequence. In this thesis we will assume throughout that i?-parity 
is conserved. 



2.3. Electro weak symmetry breaking and Higgs bosons 
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2.3 Electroweak symmetry breaking and Higgs bosons 

Electroweak symmetry breaking is caused by the fields H\ and H% acquiring vacuum expec- 
tation values 



(Hi} = 







(Ho) = 





v 2 



(2.3) 



where v± and v 2 can be chosen real and non-negative by using appropriate phases for the 
Higgs fields. They are related to the W boson mass by 



1 



m w = o5 K +^2) 



(2.4) 



and we also have the convenient expression for the Z boson mass 



2^ 



(2.5) 



where g and g' are the usual SU(2) and U(f ) gauge coupling constants. In the unitary gauge 
we then replace the fields Hi with Hi + (Hi), i = 1,2. We define the ratio of the vacuum 
expectation values, 



tan/3 = — . 

Vl 



(2.6) 



There are five physical Higgs bosons in the MSSM, H°, H 2 , H® and H^. In another 
frequently used notation, the neutral Higgs bosons are denoted by H, h and A respectively. 
We will use the notation H®, H 2 and A (and sometimes H®) for the Higgs bosons. Of the 
neutral ones, A is CP-odd and H® and H 2 are CP-even. The CP-even Higgs bosons are 
generally mixtures of the interaction eigenstates and the mixing angle is denoted by a where 
-7r/2 < a < 0. 



In Eq. (2.2) there are three parameters in the Higgs sector, mj, m 2 and B. The con- 
straints coming from minimizing the Higgs potential removes one of these and we are left 
with two independent parameters, see e.g. Ref. § for details. We have already chosen tan/3 
as one of them and it is convenient to choose the mass of the CP-odd Higgs boson, ttia, 
as our second free parameter. The masses of the other Higgs bosons are then at tree-level 
given by 



m H°H? 



m 2 H± 



(m 



,2\2 



— 4m 2 z m\ cos 2 2/3 



(2.7) 
(2.8) 



As seen by Eq. (2.7), the mass of the lightest Higgs boson m H o is at tree level bounded from 
above, 



m H o < mz I cos 2/3| . 



(2.9) 



The Higgs boson masses do however get large radiative corrections and we have used 
the renormalization group improved 2-loop leading log corrections in Ref. B. For other 
references on such effective potential approaches, see Ref. 0. The upper bound on the H® 
mass depends on the mass of the top quark, rat. For m t < 200 GeV, the upper bound is 
m H o < 150 GeV and for m t = 175 GeV, it is m H o < 130 GeV (HI. 
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2.4 Neutralinos 



The neutralinos are linear combinations of the superpartners of the gauge bosons and the 
Higgs bosons. In the basis (B, W3, H 4 , H 2 ) their mass matrix is given by 



( Mi 




V 



a Hi 
■ V2 
g V2 

n/2 






g vi 


1 S "2 


V2 


+ 


M 2 


1 gvi 

+ V2 


Sf2 

V2 


1 gvi 







gv2 

V2 








(2.10) 



The neutralino mass matrix can be diagonalized analytically to give the four neutralinos, 

(2.11) 



X° = NnB + N l2 W 3 + N l3 H° + N l4 H°, 



the lightest of which, Xi> to be called the neutralino, x, is then the candidate for the dark 
matter in the Universe. We have chosen to work with the convention where the matrix Nij 
is complex and the mass eigenvalues are all positive. The gaugino fraction of neutralino i is 
defined as 



^ = |A^ii| 2 + |^ 



/2 



(2.12) 



and we will call the lightest neutralino higgsino-like when Z g < 0.01, mixed when 0.01 < 
Z g < 0.99 and gaugino-like when Z g > 0.99 where the shorthand notation Z g = Z x g is used 
for the lightest neutralino gaugino fraction 

The neutralino mass matrix, Eq 
to dominantly quark-squark loops 



q. Q2.10 ), is valid at tree level but gets loop corrections due 
]12| , |13| ]. This effect is most important when calculating the 
relic density of higgsino-like neutralinos (as we will see in Chapter^), where both the next- 
to-lightest neutralino and the lightest chargino are close in mass to the lightest neutralino 
and in this case even small corrections to the masses are important. The most important 
one-loop corrections are corrections to entries (3, 3) and (4, 4) in the neutralino mass matrix, 
Eq. ( |2.10D , and they are given by @|3[ 

3 



'>:>,:>, 



Sn = - 



16tt 2 
3 

16tt 2 



Y b 2 m b sin(20j)Re |B (Q, b, 61) - B (Q, b, b 2 ) 
Y t 2 m t sin(2%)Re [B (Q, t, ti) - B (Q, t, t 2 )] 



where and m t are the masses of the b and t quarks, 
gm b j v gm t 



n 



\[2vnw cos f3 



and Yt — 



\/2mw sin f3 



(2.13) 
(2.14) 

(2.15) 



are the Yukawa couplings of the b and t quark, 6^ and 0^ are the mixing angles of the squark 
mass eigenstates (qi — ql cos Oq + qn sin 9q) and Bq is the two-point function for which we use 
the convention in |l^,^3). Expressions for Bq can be found in e.g. |]l4j| . For the momentum 
scale Q we use as suggested in fl2|| . Note that the loop corrections depend on the mixing 
angles of the squarks as given by the soft supersymmetry breaking parameters Ajj_d in the 
soft supersymmetry breaking potential Eq. (2.2) (see also Section 2.6 below). 



2.5 Charginos 

The chargino mass terms in the Lagrangian are given by 
W+ 



(W~H{)M i 



H 2 + 



h.c. 



(2.16) 



2.6. Squarks and sleptons 
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where the mass matrix, 

M^ = { M 1 9 J }- (2.17) 



gvi m 

is diagonalized by 

Xi = UaW- + U i2 Hi (2.18) 



,~+ 



X 



= V ll W+ + V l2 H} . (2.19) 

We choose det(C7) = 1 and U* M.y±V^ — diag(m~± , m-±) with non-negative chargino 
masses. The chargino mass matrix also gets one-loop corrections as the neutralino mass 
matrix, but these are always negligible compared to the neutralino mass corrections fl2]| and 
can hence safely be neglected. 

2.6 Squarks and sleptons 

For the squarks we choose a basis where the squarks are rotated in the same way as the 
corresponding quarks in the standard model. We follow the conventions of the Particle Data 
Group Jl5[ where the mixing is put in the left-handed <i-quark fields. 
The squark mass matrices then look like 

M 2 = ( M 2 Q + mt m u + Dt L l mt(A^^cot/3) \ 

u \ (Au- Mcot/3)m„ M 2 / + m„m^ + D RR 1 J v ' ' 

2= / K^M 2 Q K + m d ml + D d LL l m^A*, - M * tan/3) \ 
d \ (A D - Mtan/3)m d M 2 D + m*m d + ^ fi l j' 

For the sneutrinos and sleptons we in the same way get the mass matrices 

Ml = M 2 L + DZ L 1 (2.22) 

M 2 L + m e mt + D e LL l mj (A*, - fi* tan (3) 

where 



_ ivA L ' II1 e 1I1 e ' ^ LL "'eV^-B^P H ) \ (2 23) 

V (Ab -,utan/?)m e M| + mtm e + D e RR l J 



D[ L = m|cos2/3(T 3/ -e/sin 2 ^) (2.24) 
D f RR = ra|cos2/? e/ sin 2 6W (2.25) 

with T3 / being the third component of the weak isospin and e / being the charge in units of 
the elementary charge e (e > 0). In the basis we have chosen, we have 

m„ = diag(m u , m c , m t ) (2.26) 
m d = diag(TO d , m s , m b ) (2.27) 
m e = diag(m e ,m M ,m T ) (2.28) 

where we have used mt — 175 GeV for the top quark mass. 

We now need to find the mass eigenstates, fk, that diagonalize these mass matrices. The 
relations between the mass eigenstates and the interaction eigenstates Jl and fn are 

6 

ha = Y,f kT *FL' ( 2 - 29 ) 

fc=l 
6 

/*, = Y,^ T *fr (2-30) 

fc=l 
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Scan 




normal 


light 


generous 


high 


high 


light 


heavy 








Higgs 




mass 1 


mass 2 


higgsinos 


gauginos 




[GeV] 


-5000 


-5000 


-10000 


1000 


-30000 


-100 


1000 




[GcV] 


5000 


5000 


10000 


30000 


-1000 


100 


30000 




[GeV] 


-5000 


-5000 


-10000 


1000 


1000 


-1000 


1.9/i/-1.9/i 


M ma= 


' [GeV] 


5000 


5000 


10000 


30000 


30000 


1000 


2.1/i/-2.1/j 


tan /3 




1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


tan /3 




50 


50 


50 


50 


50 


2.1 


50 


m T" 


[GeV] 

























[GeV] 


1000 


150 


3000 


10000 


10000 


1000 


10000 




[GeV] 


100 


100 


100 


1000 


1000 


100 


1000 




[GcV] 


3000 


3000 


5000 


30000 


30000 


3000 


30000 


A u*n 




— 3mo 


— 3mo 


— 3mo 


— 3mo 


— 3mo 


— 3mo 


— 3mo 


^" 




3mo 


3mo 


3mo 


3mo 


3mo 


3mo 


3mo 






— 3mo 


— 3mo 


— 3mo 


— 3mo 


— 3mo 


— 3mo 


— 3mo 






3mo 


3mo 


3mo 


3mo 


3mo 


3mo 


3mo 


No. of models 


4655 


3342 


3938 


1000 


999 


177 


250 



Table 2.2: The span of MSSM parameters used for the different scans. For [i and M 2 the 
scans are uniform in the logarithms of the parameters and for the other parameters they 
are uniform in the parameters themselves. The number of models refers to the number of 
models satisfying all experimental constraints given in Section 3.1. 



where the mixing matrices T have dimension 6x3 for squarks and charged sleptons and 
dimension 3x3 for sneutrinos. 

We now have to specify the soft supersymmetry breaking parameters Ac/, Ad, Al, 
Mq, M[/, Mb, Me and Ml. To reduce the number of free parameters we make the simple 
Ansatz 

A v = diag(0,0,A t ) 

A D = diag(0, 0, A) . . 

A E = [Z - 6i) 
M Q = Mu = M D = M E = M L = m l 

which, since the matrices are diagonal, does not introduce any tree-level flavour changing 
neutral currents (FCNCs). 

2.7 GUT assumptions 

To reduce the number of free parameters further we will make the usual Grand Unified 
Theory (GUT) assumptions for the gaugino mass parameters Mi, M 2 and M3, 

Mi = - tan 2 9 W M 2 ~ 0.5M 2 (2.32) 
3 

Oi. 

M 2 = — 27r~ — M 3 ~0.3M 3 (2.33) 
sin a wct s 

where a cw is the fine-structure constant and a s is the strong coupling constant. These rela- 
tions come from the assumption that the gaugino mass parameters unify at the unification 
scale as given by the gauge coupling unification. 



2.8 Feynman rules 

The Feynman rules for the MSSM are given in Appendix The rules given there are 
basically a compilation of the rules found in Ref. [^,|^, [H| but slightly rewritten in a form 
suitable for general analytical calculations as well as numerical implementations. 



2.9. MSSM PARAMETER SCANS 
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2.9 MSSM parameter scans 

The general MSSM contains 63 free parameters |Q, but with the assumptions made in the 
previous sections in this chapter, we have reduced the number of parameters to the seven 
parameters fi, M2, tan/3, TOaj Wo, and A t . It is however a non-trivial task to sample this 
seven-dimensional parameter space in a complete way. In an attempt to do this we have 
performed several different scans in the parameter space, some of which are quite general 
and some of which are more specialized to find interesting regions of the parameter space. 



In Table 2.2 we list the different scans we have used in the explicit calculations in the 
subsequent chapters. 

Remember, though, that the actual look of our scatter plots in Chapters [| and || might 
change if different scans were used. One should especially not pay any attention to the 
density of points in different regions: it is just an artifact of our scanning. 

One might argue that the highest values of the massive parameters are unnatural and 
require fine-tuning. We have in this thesis taken a more phenomenological approach allowing 
even these high values. 
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Chapter 3 

Experimental Constraints 



For each set of parameters in the MSSM we have a unique model with given mass spectrum, 
particle properties etc. Supersymmetry is searched for both at accelerators and in dark 
matter searches and some of these models will already be excluded. In the sections below, 
both accelerator searches and dark matter searches will be discussed briefly. Note, however, 
that we only use the experimental constraints coming from accelerator searches to rule out 
models. We will, however, compare with direct dark matter searches later on. 



3.1 Accelerator searches 

Since supersymmetry introduces many new particles these can affect what is seen at accel- 
erators, either directly by finding a new particle or indirectly by changing some measured 
width or branching ratio. Below is given the currently most relevant bounds on the MSSM 
coming from accelerator searches. 



3.1.1 Neutralinos and Charginos 

The most effective limit on the chargino mass comes from LEP2, through the search for 
the process e + e~ — > x + X~ where \ + 1S the lightest chargino. This essentially puts the 
constraint m x + > \\fs (to within a few GeV p7|]). From LEP2 the present bound on the 
chargino mass is O] 

m x+ > 85 GeV. (3.1) 

The neutralinos can also be produced at LEP and would contribute to the invisible width 
of the Z boson. It is however difficult to relate this to a model-independent limit on the 
neutralino mass. In case there is no sfermion mixing, m; = mj R , and all squarks are 
degenerate in mass (except for II and in), the limit from LEP would be m x > 23 GeV when 
tan (3 > 3. We have implemented this limit by calculating the invisible width of the Z boson 
directly for each MSSM model. If the width, to which neutrinos, sneutrinos and neutralinos 
contribute, is above the experimental limit 

pmvisible < o 5Q24 Ge y (3 2 ) 

the model is excluded. 



Due to our GUT assumption, Eq. (2.32), the lightest chargino is never heavier than 



twice the neutralino mass though, and hence the chargino mass bound is more effective in 
constraining the neutralino mass in our models than the invisible Z width bound is. 
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3.1.2 Higgs bosons 

In the MSSM, the lightest Higgs boson, H®, has a tree level mass m H o < mz, but loop 



corrections can increase the mass up to about 150 GeV (as described in Section ^J). If 
no Higgs boson is seen up this mass, this would mean that the MSSM is ruled out. The 
lightest Higgs boson is searched for at LEP2, where the main processes are e + e~ — > H^Z 
and e + e~ — > H^A. The cross sections for these production channels are proportional to 
sin 2 (/3 — a) and cos 2 (/3 — a) respectively and are hence complementary to each other. More 
details about Higgs searches at LEP2 can be found in e.g. Ref. fLl| , from which we will use 



results for detection prospects of supersymmctry at LEP2 in Section 5.2 



The present LEP2 bound on the lightest Higgs boson mass is approximately given by |l8j 
m H o > 62.5 GeV, (3.3) 

but can be made more stringent by making the bound dependent on sin 2 (/3 — a). The 
bound then gets about 10 GeV higher at high sin 2 (/3 — a), but we have not included this 
more stringent mass bound here. 

3.1.3 Squarks and gluinos 

Squarks and gluinos are primarily searched for at hadron colliders. When produced they will 
eventually decay to the lightest neutralino which will escape the detector leading to a missing 
energy event. Since the squark and gluino decays depend very much on the neutralino sector 
these limits will be quite model dependent. Assuming the GUT assumptions, Eqs. ( 2.32| ) 
( p. 33 ), to hold, one can derive a limit on the squark and gluino masses as jl5| 



vrtq > 176 GeV (3.4) 
m~ g > 154 GeV. (3.5) 

There is however a controversy if there is still a window of light gluinos, ~ 1-4 GeV open 
or not. 

3.1.4 Sleptons 

Charged sleptons are searched for at e + e _ colliders where sleptons can be produced and 
eventually decay to the lightest neutralino resulting in missing energy. As for the squarks 
and gluino searches, the bounds will depend on details in the neutralino sector, especially 
on the neutralino mass. The LEP limits on the slepton masses are jl5j 

m > 37.1 GeV (3.6) 

mg > 45 GeV if m x < 41 GeV (3.7) 

m A > 45 GeV if m x < 41 GeV (3.8) 

mr > 45 GeV if m x < 38 GeV (3.9) 

3.1.5 Other searches 

Even though we have chosen the MSSM parameters to avoid tree level FCNCs, these can 
occur as one-loop corrections. It turns out that the b — > s"f decay width as measured by 
the CLEO experiment |lj]] is an important constraint on the MSSM since squark loops (in 
case of squark mixing) can change this width. We have used the following constraint on the 
decay width b — > sj, 

1.0 x 10~ 4 < BR{b -> s 7 ) < 4.0 x 10" 4 (3.10) 



3.2. Dark matter searches 
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where the branching ratio BR(b — » sj) is calculated with QCD corrections included using 
the method in Ref. 

3.2 Dark matter searches 

If neutralinos make up the dark matter in the Universe, they can also be searched for by 
different direct and indirect dark matter searches. The direct searches look for ncutralino 
scattering off nuclei in a detector. This scattering releases some energy in the detector which 
can be measured. The indirect searches look for indications of neutralino annihilation, e.g. in 
the galactic halo producing antiprotons, positrons or gamma rays or in the center of the Sun 
and Earth producing high energy neutrinos which can be detected by neutrino telescopes as 
explained in detail in Chapter |5|. 

We have not used any of these dark matter searches to exclude models, but we will 
compare the indirect detection rates in neutrino telescopes with direct detection rates in 
Chapter §. 
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Chapter 4 

Relic Density Calculations 



Since the ncutralino is a WIMP its annihilation cross section is expected to be of about 
the right magnitude to give a relic density Q, x h 2 ~ 1. The neutralino is not invented to 
solve the dark matter problem but comes from particle physics considerations and it is very 
interesting that it turns out to have a relic density in the right regime to be able to make 
up the dark matter in the Universe. 

The relic density of ncutralinos has been calculated by several authors during the years 
[fl2|, p2|-p7| and a simple, but approximate, way of calculating the relic density can be found 
in e.g. Ref. This is rather approximate since it assumes the cross section to be a nice 
function expandable in v 2 where v is the relative velocity of the annihilating particles. This 
expansion is often very bad, e.g. when there are thresholds and resonances. These problems 
have been treated in a semi-analytical way in Ref. p2[ . Instead of using these approximate 
expressions we use the full cross section and solve the Boltzmann equation numerically with 
the method given in Ref. p8[p9[] . This way we automatically take care of thresholds and 
resonances. 

When any other supersymmetric particles are close in mass to the lightest neutralino 
they will also be present at the time when the neutralino freezes out in the early Universe. 
When this happens so called coannihilations can take place between all these supersymmetric 
particles present at freeze-out. This was first noted by Griest and Seckel ||2| who investigated 
this for the rather accidental case where squarks are of about the same mass as the lightest 
neutralino. Later, coannihilations between the lightest neutralino and the lightest chargino 
were investigated by Mizuta and Yamaguchi |^5| for higgsinos lighter than the W boson. 
Drees and Nojiri |27| investigated coannihilations between the lightest and the next-to- 
lightest neutralino, which are not as important as the chargino-neutralino coannihilations. 
Recently, Drees et al. |l^] reinvestigated coannihilations for light higgsinos taking one-loop 
corrections to the neutralino and chargino masses into account. 

We have performed a more general analysis and evaluated the relic density fl x h 2 includ- 
ing coannihilation processes between all charginos and neutralinos lighter than 2.1m x for 
a general neutralino with any mass, m x , and composition, Z g . We have however not in- 
cluded coannihilations with squarks which occurs more accidentally than the in many cases 
unavoidable mass degeneracy between the lightest neutralinos and the lightest chargino. 

In the following sections, the method by which the relic density is evaluated when coanni- 
hilations are included M will be described and our results will be presented and discussed. 
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4.1 The Boltzmann equation 

We want to generalize the formulas in Ref. |28j to include coannihilations. We will do that by 
starting from the expressions in Ref. p2| which will then be rewritten into a more convenient 
form. 

Consider annihilation of N supersymmetric particles with masses m, and internal degrees 
of freedom gi. Order them such that mi < m,2 < ■■■ < t^n-X < tun- F° r the lightest 
neutralino, the notation mi and m x will be used interchangeably. The evolution of the 
number density of particle i is given by 



= -3-Hn, - J2(cT i] v tj ) (mrij - n- q n- q ) 



^[ry^-nD-T^^-nf)] (4.1) 



where 



^Xij = ^(xiX ^ X jY) (4.3) 
y 

r« = E r fc^^) (4.4) 



X 



are the total annihilation cross sections, the inclusive scattering cross sections and the in- 
clusive decay rates respectively and X and Y are (sets of) standard model particles involved 
in the interactions. The 'relative velocity' is defined by 



\/(Pi -Vjf - mf 

— (4.5) 



1 2 



EiEj 



with pi and £j being the four-momentum and energy of particle i. rij are the number 
densities of the corresponding particles given by 

d 3 Pifi (4.6) 



(2n) 



with pi being the three-momentum of particle i and fi being the equilibrium distribution 
function which in the Maxwcll-Boltzmann approximation is given by 

fi = e- E > /T (4.7) 

where T is the temperature. Since we assume that i?-parity holds, all supersymmetric 
particles will eventually decay to the LSP and we thus only have to consi der the total 
number density of supersymmetric particles n — YliLi n i- By summing Eq. (4.1) over all 
SUSY particles i we get the evolution equation for n, 

^ = -3Hn - ( n * n i - < q "j q ) ( 4 -8) 

i,3=l 



where the terms on the second and third lines in Eq. (4.1) cancel in the sum. The scattering 
rate of supersymmetric particles off particles in the thermal background is much faster than 



4.2. Thermal averaging 
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their annihilation rate, because the scattering cross sections c' Xi - are of the same order of 
magnitude as the annihilation cross sections but the background particle density nx 
is much larger than each of the supersymmetric particle densities rii when the former are 
relativistic and the latter are non-relativistic, and so suppressed by a Boltzmann factor. In 
this case, the \i distributions remain in thermal equilibrium, and in particular their ratios 
are equal to the equilibrium values, 



n n c( i 
We then get 

dn 
~dt 

where 



(4.9) 



3Hn - (a eS v) (n 2 - n 2 cq ) (4.10) 



(^ S v)=Y / {a ij v ij )^^. (4.11) 



4.2 Thermal averaging 



Now reformulate the thermal average, Eq. (4.11), into more convenient expressions. 
First, using the Maxwell-Boltzmann approximation we get |p8|,|29| 

^ = E< q = ^E^M^) ( 4 - 12 ) 

i i 

where K% is the modified Bessel function of the second kind of order 2. Then rewrite 
Eq. ( [4.1 ID as 

(a cS v) = -4" (4.13) 

where 

A -l r J W ^ 27rr2Ei{2n)32Ej - ( 4 - 14 ) 

is the total annihilation rate per unit volume at temperature T . Wij is the annihilation rate 
and is related to the cross section through)] 

= 1/;, , \ S0y = 4(jy yj (pi ■ pj) 2 - mfm 2 = 1/.', <'<j (4-15) 

where 

[s - (to, + mj) 2 ] 1/2 [s - (m* - mj) 2 ] 1/2 

Pi j = ^ . (4.16) 

For a two-body final state, Wij is given by 

W 2 r hody - 



l&iT 2 g l g j S 

lr The quantity Wjj in Ref. is Wjjji 



— J2 [ \M\ 2 dQ, (4.17) 

f internal d.o.f. 
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where k is the final center-of-mass momentum, Sf is a symmetry factor equal to 2 for 
identical final particles, and the integration is over the outgoing directions of one of the final 
particles. As usual, an average over initial internal degrees of freedom is performed. 

Now consider annihilation of two particles, i and j, with masses rrij and rrij and statistical 
degr ees o f freedom g t and gj . If we use Boltzmann statistics (good for T < to) we can put 
Eq. (4.14) into the form 



A = E/ '/<.'/.. ll <..' ' '< 



d 3 p, d 3 pj 
(2tt) 3 2E 1 (2tt) 3 2Ej 



(4.18) 



where Pi and Pj are the three-momenta and Ei and Ej are the energies of the colliding 



particles. We can now follow the procedure in Ref. []28j as done in Ref. 



some of the integrations in Eq. (4.18) to arrive at 



A = 



T 
32t? 



E 



(mi-\-mj)'- 



dsg l g j p ij W lj Ki [ — 



and perform 



(4.19) 



where K\ is the modified Bessel function of the second kind of order 1. 



Now we have what we need to perform the sum in Eq. ( 4.14 ) to get (a e gv). Let 
\ -< Pij 9i9j 



W cB = 



-Wv. 




(m, — rrij) 2 ][s — (rrii + rrijY 
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with 



Pes = pn = t:\/s - 4toi 



(4.20) 



(4.21) 



Since W jj (s) — for s < (m, + rrij) 2 , the radicand in Eq. (4.20) is never negative. 

Eq. ( 1.19] ) can be written in a form more suitable for numerical integration by using p e g 
instead of s as integration variable. From Eq. (4.21), ds = 8p e ffdp e g, and we have 



T 



We can then finally write Eq. ( 4.13| ) as 
Jo 00 dp^W.sK, 



(4.22) 



(4.23) 



This expression is very similar to the case without coannihilations, the difference being the 
denominator and the replacement of the invariant rate with the effective invariant rate. 

In the effective annihilation rate, W e ff> coannihilations appear as thresholds at -y/s equal 
to the sum of the masses of the coannihilating particles. We show an example in Fig. 4.1 
where it is clearly seen that the coannihilation thresholds appear in the effective invariant 
rate just as final state thresholds do. In Fig. 4.2 we show the differential of A with respect to 
Peff, dA/dp c ff. The Boltzmann suppression at higher p c g, contained in the exponential decay 
of Ki, is clearly visible. We have in Fig. 4.2 evaluated the modified Bessel function at the 
temperature T = m x /20 which is a typical freeze-out temperature. When the temperature 
is higher, the peak will shift to the right and when it is lower it will shift to the left. For the 
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Figure 4. 1: T he effective invariant annhiliation rate, W e g, as a function of p c s for model 1 
in Table 4.2. The final state thre shold for annihilation into W + W~ and the coannihilation 
thresholds, as given by Eq. ( 4.2C ) , are indicated. The X2X2 coannihilation threshold is too 
small to be seen. 



particular model shown in Figs. i.l-4~2|, the relic density is evaluated to be fl x h 2 = 0.030 
when coannihilations are included and il x h 2 — 0.18 when they are not. 

We end this section with a comment on the internal degrees of freedom g^. A neutralino 
is a Majorana fermion and has two internal degrees of freedom, g x o = 2. A chargino can 
be treated either as two separate species xt an d Xi > eacn with internal degrees of freedom 
g x + — g x - — 2, or, more simply, as a single species xt with g x ± = 4 internal degrees of 
freedom. 



4.3 Reformulation of the Boltzmann equation 



We can now put Eq. (4.10) into a more convenient form by instead of the number density 
considering the ratio of the number density to the entropy density, 

Y = — (4.24) 
s 

and instead of having the time t as independent variable we choose x — m\jT with m\ 
being the LSP mass and T being the temperature. By following Ref. |^8| we can write the 
evolution equation as 



/ 1/2 

dY / 7r gj mi 



dx V 45G x 2 x ' v cq 
where y cq is given by 

n on 45x 2 v-^ / ™? \ 2 „ ( n%i 



(a cS v) (Y 2 - Y 2 q ) (4.25) 



r cq = ^= Z \ T ft - KiW—\ (4-26) 



24 



Chapter 4. Relic Density Calculations 



> 



Q. 
< 




40 50 60 
Peff [C3eV] 



100 



Figure 4.2: Total differential annihilation rate p er unit volume, dA/dp c g = 
(T/TT 4 )p^ fl Ki(p c ff,T)W c s, for the same model as in Fig. 4.1. We have chosen to evaluate 
dA/dp c ff for T = m x /20 which is a typical value at freeze-out. The Boltzmann suppression 
at higher p c g should be evident. 



1/2 

with G being the Gravitational constant, the parameter gj being defined as 



.9 



1/2 



1 



T dh P 



(4.27) 



/5cff \ 3/i of f dT 

and g c ff and h c s being the effective degrees of freedom as given in the usual parameterizations 
of the energy and entropy densities 

P = .9cff(T)^T 4 , s = h cff (T)^-T 3 . (4.28) 
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We have evaluated g c s , h c ff and gV 2 using the methods given in Ref. [p8| assuming the 
QCD phase transition to occur at 150 MeV. Our results are not sensitive to the value of 
Tqcd though, since the neutralino freeze-out temperature is always much larger than Tqcd- 

To obtain the relic density we should integrate Eq. (4.25) from x — to xq = m x /To 
where To is the photon temperature of the Universe today. The relic density today in units 
of the critical density is then given by 

Q, x = p° x /pcrit = m x s Y /p clit (4.29) 

where p cri t = 3H 2 /8irG is the critical density, sq is the entropy density today and Yq is 
the solution of the integration of Eq. (4.25). With a background radiation temperature of 
T = 2.726 K we obtain 

m v 

(4.30) 



n x h z 



2.755 x W 



GeV 



4.4 Annihilation cross sections 

We have calculated all two-body final state cross sections at tree level for neutralino-neutral- 
ino, ncutralino-chargino and chargino-chargino annihilation. A complete list is given in 
Table |4T. 



4.4. Annihilation cross sections 
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Table 4.1: All two-body final states for which annihilation cross sections are calculated. The 
indices i,j,k=l, 2, 3, 4 and the indices c,d,e— 1, 2. Note that the coannihilation channels 
are only important when the mass difference is not too big. t, u and s refers to which 
channel the annihilation goes through and the particle inside the brackets is the particle in 
the propagator, u, u, d, d, v, i>, C, £, f and / are generic notations for any up- and down- 
type (s)quark, (s)neutrino, (s)lepton and (s)fermion. A sum of diagrams over (s)fermion 
generation indices and over the neutralino and chargino indices k and e is understood. 



Since we have so many different diagrams contributing, we have to use some method 
where the diagrams can be calculated efficiently. To achieve this, we classify diagrams 
according to their topology (s-, t- or u-channel) and to the spin of the particles involved. We 
then compute the helicity amplitudes for each type of diagram analytically with Reduce pi| 
using general expressions for the vertex couplings. Further details will be found in Ref. |32[ |. 

The strength of the helicity amplitude method is that the analytical calculation of a 



20 



Chapter 4. Relic Density Calculations 



given diagram only has to be performed once and the summing of the contributing diagrams 
for each given set of initial and final states can be done numerically afterwards. 



4.5 Numerical methods 

In this section we describe the numerical methods we use to evaluate the effective invariant 
rate and its thermal average, and to integrate the density evolution equation. 

We obtain the effective invariant rate numerically as follows. We generate Fortran 
routines for the helicity amplitudes of all types of diagrams automatically with Reduce, 
as explained in the previous section. We sum the Feynman diagrams numerically for each 
annihilation channel ij — > kl. We then sum the squares of the helicity amplitudes and sum 
the contributions of all annihilation channels. Explicitly, we compute 



\ - PijPki 



dW cS _ ^ 

dcos6> 2-f 2,2-nSuVs u frt. 

ijkt hclicitics 



-» kl) 



iiagrams 



(4.31) 



where 9 is the angle between particles k and i. We finally integrate numerically over cos 9 
by means of adaptive gaussian integration. 

In rare cases, we find resonances in the t- or u-channels. For the process ij — ► kl, this 
can occur when rm < and rrij > mf. at certain values of cos 9, the momentum transfer is 
time-like and matches the mass of the exchanged particle. We have regulated the divergence 
by assigning a small width of a few GeV to the neutralinos and charginos. Our results are 
not sensitive to the choice of this width, though. 

The calculation of the effective invariant rate W e & is the most time-consuming part. 
Fortunately, thanks to the remarkable feature of Eq. ( 4.23| ), W e s(p e s) does not depend on 
the temperature T, and it can be tabulated once for each model. We have to make sure 
that the maximum p e g in the table is large enough to include all important resonances, 
thresholds and coannihilation thresholds. As an extreme case, consider when the effective 
invariant rate at high p e g is 10 10 times higher than at p e g = 0. For a typical freeze-out 
temperature of T = to x /20, the Boltzmann suppression of high p e g contained in K\ in 
Eq. (4.23) results in that contributions to the thermal average from values of p e g beyond 



~ 1.5m x are negligible. For coannihilations, this value of p e ff corresponds to a mass of the 
coannihilating particle of ~ 1.8m x . To be on the safe side all over parameter space, we 
include coannihilations whenever the mass of the coannihilating particle is less than 2.1m x , 
even if typically coannihilations are important only for masses less than 1.4m x . For extra 
safety, we tabulate W c s from p e g = up to p c g = 20m x , more densely in the important low 
p e s region than elsewhere. We further add several points around resonances and thresholds, 
both explicitly and in an adaptive manner. 

To perform the thermal average in Eq. ( 4.23|) , we integrate over p e s by means of adaptive 
gaussian integration, using a spline routine to interpolate in the (p e s, W e g) table. To avoid 
numerical problems in the integration routine or in the spline routine, we split the integration 
interval at each sharp threshold. We also explicitly check for each MSSM model that the 
spline routine behaves well at thresholds and resonances. 

We finally integrate the density evolution equation ( 4.25| ) numerically from x = 2, where 
the density still tracks the equilibrium density, to xq — m x /T^. We use an implicit trape- 
zoidal method with adaptive stepsize. The relic density at present is then evaluated with 
Eq. (|O0l) . 

A more detailed description of the numerical methods will be found in a future publication 

M. 



4.6. Results 
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Figure 4.3: Neutralino relic density including neutralino and chargino coannihilations 
versus a) neutralino mass m x and b) neutralino composition Z g /(1 — Zg). The horizontal 
lines indicate the cosmologically interesting region 0.025 < Q x h 2 < 1. 




Figure 4.4: Ratio of the neutralino relic densities with and without neutralino and chargino 
coannihilations versus a) neutralino mass m x and b) neutralino composition Z g /{1 — Zg). 



4.6 Results 



We now present the results of our relic density calculations for all the models in Table 2.2. 
We will focus on the effect of coannihilations, since this is the first time they are included 
for general neutralino masses and compositions. 

Fundamentally, we are interested in how the inclusion of coannihilations modifies the 
cosmologically interesting region and the cosmological bounds on the neutralino mass. We 
define the cosmologically interesting region as 0.025 < tt x h 2 < 1. In this range of Sl x h 2 the 
neutralino can constitute most of the dark matter in galaxies and the age of the Universe is 
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Table 4.2: Some representative models for which coannihilations are important (examples 
1-5) and one model (example 6) for which they are not. We give the seven model parameters, 
the masses of the lightest neutralinos and of the lightest chargino, the gaugino fraction of 
the lightest neutralino and the relic densities without coannihilations included and with. 



long enough to be compatible with observations (see Chapter [I]). The lower bound of 0.025 
is somewhat arbitrary, and even if Q, x h 2 would be less than 0.025 the neutralinos would still 
be relic particles, but only a minor fraction of the dark matter in the Universe. 

We start with a short general discussion and then present more details in the following 
subsections. 

Fig. 4.3 shows the neutralino relic density Q, x h 2 with coannihilations included versus the 
neutralino mass m x and the neutralino composition Z g /(1 — Z g ), respectively. The lower 
edge on neutralino masses comes essentially from the LEP bound on the chargino mass, 
Eq. ©. 

The neutralino is a good dark matter candidate in the cosmologically interesting region 
limited by the two horizontal lines. There are clearly models with cosmologically interesting 
relic densities for a wide range of neutralino masses and compositions. The cosmologically 



interesting region will be discussed more in Section 4.6.5 



Coannihilations increase the 
< (ciiivii). Recalling that 



The effect of neutralino and chargino coannihilations on the value of the relic density 
is summarized in Fig. 4.4, where we plot the ratio of the neutralino relic densities with 
and without coannihilations versus the neutralino mass m x and the neutralino composition 
Z g /(1 — Zg). In many models, coannihilations reduce the relic density by more than a factor 
of ten, and in some others they increase it by a small factor, 
relic density if the effective annihilation cross section (er ffv) 
(a e gv) is the average of the coannihilation cross sections (see Eq. (4.11)), this occurs when 
most of the coannihilation cross sections are smaller than (anvn) and the mass differences 
of the coannihilating particles are small. 

Table 4.2 lists some representative models where coannihilations are important plus 
one model where coannihilations are negligible. Example 1 contains a light higgsino-likc 
neutralino, example 2 a heavy higgsino-like neutralino. Examples 3 and 4 have ~ \Mi\, 
and example 5 has a very pure gaugino-like neutralino. Example 6 is a model with a 
gaugino- like neutralino for which coannihilations are not important. 

In Fig. 4.5 we show the reduction in relic density due to the inclusion of coannihilations as 
a function of \fj,/Mi \ . A rule of thumb is that coannihilations are important when \fi/Mi | < 2. 
But exceptions are found, as can be seen in Fig. 4.5. Notice that when \fj,/M±\ <C 1, the 
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Figure 4.5: Ratio of the relic densities with and without coannihilations versus \n/M\\. 
Coannihilations are important when \n/M\\ < 2. 



neutralino is higgsino-like; when |/x/Mi| ^> 1, the neutralino is gaugino-like; and when 
\li/M\\ ~ 1, the neutralino can be higgsino-like, gaugino-like or mixed. 

It can be convenient to have a criterium for when coannihilations are important in terms 
of the composition as well. A rule of thumb is that coannihilations are important when 
Z g < 0.23 for m x < 200 GeV and when Z g /(1 ~ Z g ) < (m x /300 GeV) 3 for m x > 200 GeV. 
There are exceptions to this rule as well, as can be seen in Fig. [O^ where the ratio of relic 
densities with and without coannihilations is plotted versus the neutralino mass, the left 
panel for points satisfying the present criterion, the right panel for those not satisfying it. 

In the following subsections, we present the cases where we found that coannihilations 
are important and explain why. We first discuss the already known case of light higgsino- 
like neutralinos, continue with heavier higgsino-like neutralinos, the case ~ \Mi\ and 
finally very pure gaugino-like neutralinos. We then end this section by a discussion of the 
cosmologically interesting region. 



4.6.1 Light higgsino-like neutralinos 

We first discuss light higgsino-like neutralinos, m x < mw, Z g < 0.01, since coannihilation 
processes for these have been investigated earlier by other authors [^2|, p5|, p7fl . 

Mizuta and Yamaguchi |2j| stressed the great importance of including coannihilations 
for higgsinos lighter than the W boson. For these light higgsinos, neutralino-neutralino 
annihilation into fermions is strongly suppressed whereas chargino-neutralino and chargino- 
chargino annihilations into fermions are not. Since the masses of the lightest neutralino 
and the lightest chargino are of the same order, the relic density is greatly reduced when 
coannihilations are included. Mizuta and Yamaguchi claim that because of this reduction 
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Figure 4.6: Ratio of the relic densities with and without coannihilations versus neutralino 
mass m x . Coannihilations are generally not important when Z g > f(m x ), with / given in 
the text. 



light higgsinos are cosmologically of no interest. 

Drees and Nojiri |27j included coannihilations between the lightest and next-to-lightest 
neutralino, but did not include those between the lightest neutralino and chargino, which 
are always more important. In spite of this, they concluded that the relic density of a 
higgsino-like neutralino will always be uninterestingly small unless m x > 500 GeV or so. 

Drees at al. JlJ] then re-investigated the relic density of light higgsino-like neutralinos. 
They found that light higgsinos could have relic densities as high as 0.2, and so be cosmo- 
logically interesting, provided one-loop corrections to the neutralino masses are included. 

We agree with these papers qualitatively, but we reach different conclusions. We show 
our results in Fig. 4.7, where we plot the relic density of light higgsino-like neutralinos versus 
their mass with coannihilations included, as well as the ratio between the relic densities with 
and without coannihilations. The Mizuta and Yamaguchi reduction can be seen in Fig. 4.7b 
below 100 GeV, but due to the recent LEP2 bound on the chargino mass the effect is not 
as dramatic as it was for them. If for the sake of comparison we relax the LEP2 bound, the 
reduction continues down to 10~ 5 at lower higgsino masses and we confirm qualitatively the 
Mizuta and Yamaguchi conclusion — coannihilations are very important for light higgsinos 
— but we differ from them quantitatively since we find models in which light higgsinos have 
a cosmologically interesting relic density. For the specific light higgsino models in Drees 
et al. we agree on the relic density to within 20-30%. We find however other light 
higgsino-like models with higher fl x h 2 ~ 0.3, even without including the loop corrections to 
the neutralino masses. 



4.6. Results 



31 




Figure 4.7: For higgsino-likc ncutralinos (Z g < 0.01) we show a) the relic density with coan- 
nihilations included and b) the ratio of the relic densities with and without coannihilations 
versus the neutralino mass. The horizontal lines in a) limit the cosmologically interesting 
region 0.025 < VL x h 2 < 1. 



So there is a window of light higgsino models, m x ~ 75 GeV, that are cosmologically 
interesting. All these models have tan/3 < 1.6 and those with the highest relic densities have 
tan/3 ~ 1.2. These models escape the LEP2 bound on the chargino mass, m x + ~ 85 GeV, 
because for tan (3 < 2 the mass of the lightest neutralino can be lower than the mass of the 
lightest chargino by tens of GeV. By the same token, coannihilation processes are not so 
important and the relic density in these models remains cosmologically interesting. Most of 
these models will be probed in the near future when LEP2 runs at higher energies, but some 
have too large a chargino mass (m+ > 95 GeV) and too large an boson mass (m H o > 90 
GeV) to be tested at LEP2. Thus ~ 75 GeV higgsinos with tan/3 < 2 may remain good 
dark matter candidates even after LEP2. 



4.6.2 Heavy higgsino-like neutralinos 

Coannihilations for higgsino-like neutralinos heavier than the W boson have been mentioned 
by Drees and Nojiri j27), who argued that they should not change the relic density by much, 
and by McDonald, Olive and Srednicki |2J], who warn that they might change it by an 
estimated factor of 2. We typically find a decrease by factors of 2-5, and in some models 



even by a factor of 10 (see the right hand side Fig. 4.7b). 

For m x > m\y, the lightest and next-to-lightest neutralinos and the lightest chargino are 
close in mass, and they annihilate into W bosons besides fermion pairs. While the annihi- 
lation and coannihilation cross sections into W pairs are comparable, the coannihilation of 
X1X2J X1X1 an d X2X1" m t° fermion pairs is stronger than the X1X1 ~* ff annihilation which 
is suppressed. This gives the increase in the effective annihilation rate that we observe. 

As a result, the smallest and highest masses for which higgsino-like neutralinos heavier 
than the W boson are good dark matter candidates shift up from 300 to 450 GeV and from 
3 to 7 TeV respectively. 

Together with the result in the previous subsection, we conclude that higgsino-like neu- 
tralinos (Z g < 0.01) can be good dark matter candidates for masses in the ranges 60-85 
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GeV and 450-7000 GeV. 



4.6.3 Models with ~ |Mi| 

Coannihilations for mixed or gaugino-like neutralinos have not been included in earlier cal- 
culations. It has been believed that they are not very important in these cases. On the 
contrary, when ~ \M%\ and m x > mw there is a very pronounced mass degeneracy 
among the three lightest neutralinos and the lightest chargino. The ensuing coannihilations 
can decrease the relic density by up to two orders of magnitude or even increase it by up to 



a factor of 3. This is easily seen in Fig. 4.5 as the vertical strip at \[i/Mi \ ~ 1. 

If the lightest neutralino is mixed, Z g ~ 0.5, coannihilations can increase the relic density, 
whereas if it is more higgsino-like or gaugino-like they will decrease it. This because the 
annihilation cross section for mixed neutralinos is generally higher than those for higgsino- 
like or gaugino-like neutralinos. 

The largest decrease we see for this kind of models is when \M\\ is slightly less than 
and both are in the TeV region. In this case, the lightest neutralino is a very pure bino, and 
its annihilation cross section is very suppressed since it couples neither to the Z nor to the 
W boson. The chargino and other neutralinos close in mass have much higher annihilation 
cross sections, and thus coannihilations between them greatly reduce the relic density. This 
big reduction suffices to lower Q, x h 2 to cosmologically acceptable levels if Z g < 0.96. This 
reduction does not occur for masses much lower than a TeV, because the terms in the 
neutralino mass matrix proportional to the W mass prevent such pure bino states and the 
severe mass degeneracy. 

To conclude, when ~ |Afi|, coannihilations are very important no matter if the 
neutralino is higgsino-like, mixed or gaugino-like. The relic density can be cosmologically 
interesting for these models as long as the gaugino fraction Z g < 0.96: these neutralinos are 
good dark matter candidates. 



4.6.4 Gaugino-like neutralinos with \n\ ^> \M±\ 

When ^> \M\ \, the lightest neutralino is a very pure gaugino. According to the GUT 
relation Eq. ( 2.32| ), the supersymmetric particles next in mass, the next-to- lightest neutralino 



and the lightest chargino, are twice as heavy. So we expect that coannihilations between 
them are of no importance.^] In fact, as discussed in section 4.5, coannihilations would need 



to increase the effective cross section by several orders of magnitude for these large mass 
differences. 

This actually happens in some cases. They show up as the small spread at high |/i/Mi| 
in Fig. 4.5. In these models, the lightest neutralino is a very pure bino (Z g > 0.999) and the 



squarks are heavy. Its annihilation to fermions is suppressed by the heavy squark masses, 
and its annihilation to Z and W bosons is either kincmatically forbidden or extremely 
suppressed because a pure bino does not couple to Z and W bosons. On the other hand, 
the lightest chargino annihilates to gauge bosons and fermions very efficiently. The huge 
increase in the effective cross section, compensated by the large mass difference, reduces the 
relic density by 10-20%. However, the relic density before introducing coannihilations was 
of the order of 10 3 -10 4 , and this small reduction is not enough to make these special cases 
cosmologically interesting. 
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Figure 4.8: Neutralino masses m x and compositions Z g /(1 — Z g ) for cosmologically inter- 
esting models a) with and b) without inclusion of coannihilations. 

4.6.5 Cosmologically interesting region 



We now summarize when the neutralino is a good dark matter candidate. Fig. 4.8 shows 
the cosmologically interesting region 0.025 < fl x h 2 < 1 in the neutralino mass-composition 
plane Z g /(1 — Z g ) versus m x . 

The light higgsino-like region does not extend to the left and down due to the LEP2 
bound on the chargino mass. The lower edge in gaugino fraction at Z g ~ 10~ 5 is the border 
of our survey (how high \Mi\ is allowed to be). The upper limit on Z g and the upper limit 
on the neutralino mass come from the requirement fl x h 2 < 1. The hole for higgsino-like 
neutralinos with masses 85-450 GeV comes from the requirement fi x h 2 > 0.025. 

We see that coannihilations change the cosmologically interesting region in the following 
aspects: the region of light higgsino-like neutralinos is slightly reduced and the big region 
of heavier higgsinos is shifted to higher masses, the lower boundary shifting from 300 GeV 
to 450 GeV and the upper boundary from 3 TeV to 7 TeV. 

The fuzzy edge at the highest masses is due to models in which the squarks are close in 
mass to the lightest neutralino, in which case t- and it-channel squark exchange enhances 
the annihilation cross section. In these rather accidental cases, coannihilations with squarks 
are expected to be important and enhance the effective cross section even further. Thus, 
the upper bound on the neutralino mass of 7 TeV is an underestimate. 



2 In models with non-universal gaugino masses, the lightest gaugino-like neutralino can be almost degen- 
erate with the lightest chargino, and coannihilations can be important, as examined e.g. in Ref. |33| 
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Chapter 5 

Neutralino Detection by 
Neutrino Telescopes 



As mentioned in Chapter [l], neutralinos can accumulate in the Sun p^-p6[| or the Earth p7| ; 
annihilate and produce detectable muon neutrinos. In this chapter, the method we have used 
to predict the muon fluxes resulting from neutralino annihilation in the Sun or Earth will 
be described. We will also discuss, what signal flux levels that can be probed by neutrino 
telescopes and what a detected signal will tell us about the neutralino mass. 

There are many neutrino telescopes in use, e.g. Baksan |3q|, MACRO [p9f, Kamiokandc 
[ ^0| , and others being built and/or proposed, Amanda |42||, Nestor |||, and maybe 
others. A neutrino telescope consists of water or ice situated well below the ground (to 
minimize the background coming from atmospheric muons). When a neutrino passes by it 
may interact with the ice or rock surrounding the detector and produce a lepton. If the 
neutrino is a muon neutrino and hence the lepton a muon, it will neither decay too fast 
or get stopped too fast and may travel several kilometers (depending on its energy) before 
getting stopped. As the muon moves through the water (ice) it will emit Cerenkov radiation, 
which can be detected by photomultipliers. From this signal, the direction of the muon and 



thus the muon neutrino can be reconstructed. In Fig. 5.1 we show the principle of a neutrino 
telescope as described above. 

In the following sections, the steps performed to obtain a prediction of the muon flux for 
a given set of MSSM parameters are described. Note however that except for Section 5.3, 
the results obtained in this chapter are valid for any WIMP and not only neutralinos. 



5.1 Neutralino capture and annihilation rates 
5.1.1 Capture and annihilation rates 

When the Sun and Earth moves through our galactic halo, neutralinos may scatter off 
nuclei in them and lose enough energy to get gravitationally trapped |34j. They will then 
oscillate back and forth, occasionally scatter and after a while accumulate in the center of 
the Sun |35],|3(| and the Earth |37|] where they can annihilate. The evolution equation for 
the number of neutralinos, N, in the Sun or the Earth is given by 

^- = C - C A N 2 ~ C E N (5.1) 

where the first term is the neutralino capture, the second term is twice the annihilation rate 
Ta = \CaN 2 and the last term is neutralino evaporation. The evaporation term can be 
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Figure 5.1: Schematic view of a neutrino telescope like Amanda. The dashed lines represent 
the Cerenkov cone of light emitted by the muon when it traverses the ice. The Cerenkov 
light is picked up by the photomultipliers from which the muon track can be reconstructed. 
Note that the muon scattering angle is exaggerated in the figure. 



neglected for neutralinos heavier than about 5 GeV Q,|4|] and since we are not interested 



in these low- mass neutralinos we can safely drop the last term in Eq. (5.1). If we solve 



Eq. (5.1) for the annihilation rate Ta we get 



1~. ,2* 



T A = -Ctanh 2 -, r = yJCC A (5.2) 
2 T 

where r is the time scale for capture and annihilation equilibrium to occur. In most cases 
where the muon fluxes are within reach of present and near-future telescopes, equilibrium 
will have occurred and the annihilation rate is at 'full strength', Fa — \C. Note that in this 
case, the annihilation rate is determined by the elastic scattering cross sections, on which 
C depends, and not by the annihilation cross section. In our calculation we have of course 
used the full expressions without assuming that the annihilation occurs at 'full strength'. 

The capture rate C depends on, among other things, the local halo mass density, p Xl the 
velocity dispersion of dark matter particles in the halo, v — \J (v 2 ), the elastic scattering 
cross sections and the composition of the Earth and Sun and we have used the convenient 
expressions given in Ref. |i] based on the formulas in Ref. f46| . The main uncertainties in 
the capture rate come from the local halo mass density, p x , which is uncertain of about a 
factor of two or so, and the velocity dispersion v. We have chosen p x — 0.3 GeV/cm 3 and 
v = 270 km s _1 . Estimates of p x can be found in e.g. Ref. Jj^ . 

Note that if one-loop corrections to the neutralino coupling to Higgs bosons Q are 
included, which we have not, this coupling can for higgsino-like neutralinos either increase 
by more than two orders of magnitude or in accidental cases be reduced to exactly 0. 
This means that the spin-independent scattering cross sections for higgsinos can get greatly 
increased or reduced which mainly effects the capture rate in the Earth (and direct detection 
experiments which we don't discuss here). For mixed or gaugino-likc neutralinos these one- 
loop corrections are expected to be small. 



5.1. Neutralino capture and annihilation rates 



37 




Figure 5.2: Projected angular distributions of neutralino generated neutrinos from the Earth 
for neutralino masses of 50 GeV (solid), 100 GeV (dashed), 200 GeV (dotted), 350 GeV 
(dash-dotted) and 750 GeV (wide dotted). The analogous distributions from the Sun are 
simply narrow peaks at 9 V = 0° 



5.1.2 Annihilation profiles 

The annihilation rate per volume clement is given by 

T = n 2 (aAv) 



(5.3) 



where n x is the number density of neutralinos and {<jav) is the thermally averaged annihi- 
lation cross section. The number density of neutralinos is given by |4q,[49| 



n(r) = n(0)e- r2/2 < 



with 



3fcT 



-l 1/2 



0.56i? s 



(5.4) 



(5.5) 



where we in the last step have assumed that the annihilation takes place in the Earth. We 
have used the Earth radius i?® = 6378 km, the core temperature T ~ 6000 K and the core 
density p ~ 13.1 g cm~ 3 . In Fig. 5.2 the resulting projected neutrino angle is shown for 
neutralino annihilation in the Earth. 

For the Sun, the annihilation region is very concentrated to the core and subtends a 
negligible solid angle as seen from a neutrino telescope at the Earth. 
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5.2 Muon fluxes - Monte Carlo simulations 

A detailed Monte Carlo simulation of the muon fluxes at neutrino telescopes for different 
annihilation channels and neutralino masses was first done by Ritz and Seckel fl50| ]. We 
have followed their approach, except for neutrino propagation through the Sun where we 
have used detailed Monte Carlo simulations ffl,E4| instead of the approximate analytical 



formulae in Ref. |30|. With respect to calculations using Ref. |50| (e.g. Ref. |52|), this Monte 
Carlo treatment of the neutrino propagation through the Sun increase the muon fluxes by 
5-20%. 

5.2.1 Annihilation channels and branching ratios 

From the previous section we know how to calculate the neutralino annihilation rates and 



it is straightforward (by e.g. the same methods as described in Section 4.4) to calculate the 
annihilation branching ratios into different annihilation channels. Since the temperature in 
both the center of the Earth and Sun is so small, the neutralinos are highly non-relativistic 
and we can to a good approximation use the zero relative velocity annihilation cross sections 
to calculate the branching ratios. The annihilation channels of any significance for the muon 
fluxes are cc, bb, ti, t+t~, W+W~, Z°Z°, Z°H°, ZH$, H?H§, H$H° and H ± W ± . Note 
that the branching ratio into neutrinos directly is zero in the non-relativistic limit and hence 
the only neutrinos we get are those coming from decay of other annihilation products. We 
can also not have annihilation into Z°H^ or 2 H® 2 in this limit since the initial state is CP- 
odd and so must the final state be. Lighter quarks will not contribute since the annihilation 
cross section into fermions is approximately proportional to the mass of the fermion squared. 
Of the charged leptons, only muons and tauons are interesting as potential muon neutrino 
producers but as we will see in the next subsection, muons will be stopped well before they 
decay whereas tau leptons do decay before getting stopped so the only lepton channel we 
have to consider is t + t~. 

5.2.2 Charged lepton interactions in the Earth 

For relativistic charged particles (other than electrons) the mean energy loss is given by the 
Bcthc-Bloch equation |l5[ 



dE „ 2 Z 1 



Kz 



dx A(3 2 
where 



1 2m e c 2 /3 2 7 2 T max 2 hw p 1 
-In — ft / -hi/?7+ - 



(5.6) 



K = A-KN A rlm e c 2 ~ 0.307 MeV cm 2 mob 1 (5.7) 

2m e c 2 /3 2 7 2 

Tmax = 4 + 2 7 m e /Af + (m e /Mf (5 ' 8) 

huj p = 28M6y/p(Z/A) eV (5.9) 

with M being the mass of the particle, fie its velocity, Z the atomic number of the stopping 
material, A its atomic weight, p its density and / the mean excitation energy which is given 
approximately by 

7~16Z a9 eV. (5.40) 



At high energies (> 500 GeV for muons in ice), Eq. (5.6) is only a lower limit for the energy 
loss due to that radiative effects not included in Eq. (5.6) start dominating at high energies. 
This is however not important in our present analysis, but will be included in Section 5.2.9 
where more accurate expressions are needed. 



5.2. Muon fluxes - Monte Carlo simulations 
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We want to see if muons and tau leptons have time to decay before they get stopped. For 
muons in the core of the Earth (which consists of mainly liquid iron of density 13.1 g/cm 3 ) 



the minimum of Eq. (5.6) is 

Earth 



dE\ 
~dt) 



5.7xl0 8 GeV/s (5.11) 



from which the mean stopping time is given by r s to P = E/(dE/dt). Since the decay time is 
given by 7Td 0C = Et^/M we find that when the ratio 

= (5.12) 

7Tdcc § r doc 

is much smaller than one, the particles will get stopped well before they have time to decay. 
For muons we get 

) Earth 
<8.4xl0" 5 (5.13) 

and hence any muon produced will be stopped before having time to decay producing any 
muon neutrinos. For tau leptons, on the other hand, we can in the same way find that the 
maximal energy loss per second (up to E T — 5000 GeV) is given by 



dE\ Earth 

J max 



9.3 x 10 s GeV/s (5.14) 



from which it is seen that the upper limit on the energy loss of a tau lepton, ^T^dE / dt is 
less than a GeV even for TeV energy tauons. Hence energy loss of tau leptons can safely be 
neglected. 

Note that the simple estimates in this subsection could change by a factor of 2-3 or so 
at high energies, > 500 GeV, if radiative effects are included. The conclusion that muons 
get stopped and that tau lepton interactions can be neglected would be the same though. 



5.2.3 Charged lepton interactions in the Sun 



For charged lepton interactions in the core of the Sun, which is a plasma, Eq. (5.6) has to 
be replaced by 153] 



dE 
d.v 



4^eo {Pcf 



■In 



Am e c 2 7/3 2 



(5.15) 



where lo p is the plasma frequency as given by Eq. (5.9) and A is a number of order unity 
which we put equal to 1. By using that the composition of the core of the Sun is 24% 1 H and 
64% 4 He with a density of 148 g/cm 3 |54j we then get the minimal energy loss for muons in 
the core of the Sun to be 



-^L ) ~ 7.5 x 10 9 GeV/s 
dt ' 



Sun 



and hence 

/ 7"stop 
V 7 r dcc 



Sun 



< 6.4 x 10~ 6 . 



(5.16) 



(5.17) 
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Muons are thus stopped well before they decay and can hence be considered as absorbed. 
For tau leptons we get the maximal energy loss (up to E T = 5000 GeV) to be 

Sun 

~ 1.4 x 10 10 GeV/s (5.18) 

max 

which implies that the energy loss is only a few GeV even for TeV tau leptons. Hence tau 
lepton energy loss can be neglected. 

To conclude the previous subsection and this one, muons get stopped well before they 
decay in both the Earth and Sun and the energy loss of tau leptons can be neglected in both 
the Earth and Sun. 




5.2.4 Heavy hadron interactions in the Sun and Earth 

The heavy quarks produced in neutralino annihilation will form mesons and baryons which 
may interact before they decay. The top quarks will decay before they even have time to 
form any hadrons and their interactions with the surrounding medium can be neglected. 
For c and b quarks however, interactions with the surrounding medium has to be taken care 
of. This is done in an approximate fashion where the decay/hadronization is simulated as if 
in vacuum as described below and the interactions that may have occurred are introduced 
afterwards as a general energy decrease. This approximation is reasonable when the number 
of interactions is not more than a few which is the case for moderately heavy neutralinos, 
m x < 500 GeV. For heavier neutralinos neither the cc nor the bb channel will dominate and 
hence the approximation is justified. 

The cross sections for c and b hadron scattering off a nucleon can be estimated by noting 
that the scattering cross section for any hadron off a proton is approximately given by [pq] 

(5.19) 

where (r 2 t )hadron is the mean squared strong interaction radius of the hadron. Povh et al. 
found that (r s 2 t ) A ~ 0.58±0.02 fm 2 , (r 2 ^ ~ 0.41 ±0.02 fm 2 and (r 2 t )j/* ~ 0.04±0.02 fin 2 . 
If we assume that the decrease in mean squared radius is constant for each light quark we 
change to a c or b quark (justified by experiments) we find that the mean squared strong 
interaction radii for c or b mesons and baryons are given by 

(r s 2 t)c/6- mC so„ ^ 0.23 fm 2 (5.20) 

(r s 2 t > c / 6 -baryon - 0.40 fm 2 (5.21) 



which together with Eq. ( |5.19| ) yields 

CTc/fc-mcson ^ 14 mb (5.22) 
°"c/h-baryon ~ 24 mb ( 5 - 23 ) 

From the simulations it is known how far the leading hadron has moved before decaying 
and the probability that it should have undergone one or more interactions is then easily 
calculated given the cross sections above. If an interaction is found to have occurred the 
leading hadron of the new jet takes the energy fraction z of the initial hadron. The average 
energy transfers, (z), used are those calculated by Ritz and Seckel IpCfl , 

(z) ~ 0.6 — - for c hadrons (5-24) 

mi 

(z) ~ 0.7 for b hadrons (5.25) 

where m, is the mass of the initial c hadron and m c is the mass of the c quark. The same 
amount of energy decrease is assumed to apply to the produced neutrino. For the bb channel, 
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interactions are only significant in the Sun and only for heavier neutralinos, m x > 50 GeV. 
The effect is very dramatic for even more massive neutralinos but for those other annihilation 
channels will dominate. 



5.2.5 Monte Carlo simulations 

Of the annihilation channels mentioned above, gauge bosons and tau leptons can decay di- 
rectly to neutrinos but the quarks will hadronize and eventually give rise to muon neutrinos. 
The Higgs bosons will decay mainly to quarks which will hadronize as well. All decays and 
hadronizations are simulated with the Lund Monte Carlo Jetset 7.4 and Pythia 5.7 Q 
for each of the annihilation channels cc, bb, ti, t + t~, W + W~ and Z°Z° for the different 
neutralino masses m x = 10, 25, 50, 80.2, 91.3, 100, 150, 175, 200, 250, 350, 500, 750, 1000, 
3000 and 5000 GeV. Note that the annihilation channels containing Higgs bosons do not 
need to be simulated separately since the Higgs bosons decay to particles contained in these 
six 'fundamental' channels mentioned above and their contribution to the muon neutrino 
flux can thus be calculated as soon as the Higgs masses and their decay channels are known. 
For each mass and annihilation channel, 2.5 x 10 5 events have been simulated and all muon 
neutrinos produced have been kept. Hence, a neutrino flux is obtained for any of the given 
annihilation channels and neutralino masses. 



5.2.6 Neutrino interactions and cross sections 

The neutrino-nucleon charged and neutral current cross sections are approximately given by 

ctcc - aE v (5.26) 
a NC ~ bE v (5.27) 

where the coefficients a and b are given by 



(5.28) 
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(5.29) 



where GRV structure functions |p7j have been used down to Q 2 — 0.3 GeV 2 . These cross 
sections agree well with neutrino experiments on isoscalar targets as well as with those 
obtained using other structure functions, like CTEQ2D |58|]. 



5.2.7 Neutrino interactions in the Sun 

Ritz and Seckel [50| ] considered neutrino interactions on the way out of the Sun in an approx- 
imate way where neutral current neutrino-nucleon interactions were assumed to be much 
weaker than charged current interactions and the energy loss was assumed to be continuous. 
Neither of these approximations are very good and hence we have instead simulated the 
neutrino interactions on the way out of the Sun with Pythia 5.7 ^6|. For the Earth, the 
effective thickness is not big enough to be of any importance and neutrino energy loss or 
absorption on the way to the detector can thus be neglected. 
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Figure 5.3: The neutrino spectrum (v^ + v^) weighted by E„ at the surface of the Sun for 
the W + W~ channel with a neutralino mass of 500 GeV. The two peaks in the Ritz and 
Seckel spectrum correspond to and respectively. 



The effective thickness for the Sun is calculated by using the solar model in Ref. ]54[ |. 
The effective thickness of the Sun is given by RqIq of which the part in hydrogen is RqIh 
where i?© is the radius of the Sun. The J's are given by 

I = [ p(x)dx ~ 21.4 g cm" 3 (5.30) 
Jo 

I H = [ p{x)X H (x)dx~ll.0gcm.- 3 (5.31) 
Jo 

where x = r/i?©, p{x) is the solar density and Xh{x) is the hydrogen mass fraction. In 
terms of protons and neutrons the corresponding integrals would be 

I P = \{Io + Ih) = 16.2 gem- 3 (5.32) 

/„ = i(/ -/ff) = 5.2gcm- 3 . (5.33) 



With these effective thicknesses and the neutrino-nucleon cross sections, Eqs. (5.26) 



( 5.29 ), at hand it is then straightforward to calculate the probability that a given (anti)neu- 
trino has participated in an interaction (charged or neutral current), and if it has (and the 
interaction is a neutral current interaction) simulate it with Pythia and proceed with the 
same procedure until the neutrino has reached the surface of the Sun. 

In principle one should also take into account the fact that all annihilations don't occur 
exactly at the center of the Sun which will introduce a smearing of the effective thickness. 
This effect can be estimated jl8| and is found to be very small. Hence it is a very good 
approximation to assume that all neutrinos from the Sun originate from the center. 

The difference between this approach and the Ritz and Seckel approach is shown in 



Fig. 5.3 where the neutrino flux weighted by the neutrino energy squared is shown. Note 
that E^dN v jdE v is approximately proportional to the muon flux since both the neutrino- 
nucleon cross sections and the muon range are approximately proportional to the energy 



5.2. Muon fluxes - Monte Carlo simulations 



43 



of the neutrino/muon. The mean energy of the neutrinos at the surface of the Sun is 
about the same with the Ritz and Scckcl approach and this more detailed analysis but the 
distribution is different and since the muon flux is proportional to the second moment of the 
distribution (as explained above) one should expect a difference in the predicted muon fluxes 
at a detector. In fact, the total muon flux with this method is about 5-20% higher than 
with the Ritz and Seckel approach (with the higher difference at higher masses). Except 
for this difference in neutrino interactions our results agree well with the Ritz and Seckel 
results |5^| as well as with the analytical results in Ref. [ |59"| ]. 



5.2.8 Neutrino interactions at the detector 

When a neutrino comes close to the detector it may interact and produce a muon via a 
charged current interaction. This process is also simulated with Pythia 5.7 ]5(| where 
information of not only the muon energy but also the muon angle with respect to the 
neutrino is kept. For neutralino annihilations in the Earth, the size of the annihilation 
region has also been included according to the distributions in Eqs. (5. 



,-(|l). 



5.2.9 Muon interactions 



When a muon is produced it can travel several kilometers (depending on energy) before 
reaching a detector. When the muons travel through the ice or rock surrounding the detector 
they may interact and lose energy, where the energy loss is approximately given by 



dx 



-a - j3E^ 



where the coefficients a and /3 are fitted to the energy losses calculated in Ref. 
given by 
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for muons propagating in water or ice and 
(a/p) Iock ~ 0.00221 GeV g- 1 
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cm 
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(5.34) 



and are 



(5.35) 



(5.36) 



for muons propagating in rock. The errors of these parameterizations are less than 2% in the 
region 30-10000 GeV and less than 6% in t he re gion 10-30 GeV. Hence they are sufficiently 
accurate for our needs. By integrating Eq. (5.34) we get the mean energy of the muons after 
having traversed a distance x of the detector surroundings to be 



EM 



-fix 



(5.37) 



where E® is the initial muon energy. From this relation we find that the range of a muon of 



energy E° is 



a/0 



(5.38) 



which for low energies, E® < a/(3 ~ 500 GeV, is approximated by = E®/a. Note that 
the radiative effects contained in the /3-term in Eq. (5.34) are actually stochastic and will 
lead to energy straggling, i.e. some muons will lose more energy and some less. This is 
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however only important when > a/ ft ~ 500 GeV and for our needs the use of the mean 



energy as given by Eq. ( 5.37 ) is good enough. 

The muons will also undergo multiple Coulomb scattering on their way to the detector in 
which process they don't lose energy but the angular distribution of the muons gets smeared. 
For the multiple Coulomb scattering we have used the formulas in Ref. jl5) . 

Since the neutrino flux to a very good approximation is constant in the region surrounding 
the detector where neutrino-interactions producing detectable muons occur, we for each 
produced muon choose a distance between and L M away from the detector where the muon 



was produced and degrade its energy on its way to the detector according to Eq. (5.37). We 
also take care of the Multiple Coulomb scattering occurring during this passage of matter 
as described above. 

5.2.10 Resulting muon fluxes 

With the methods described above, we have what we need to calculate the muon flux at a 
detector for a given neutralino (or any WIMP) mass and annihilation channel. As described 
above, 2.5 x 10 5 annihilations are simulated for each neutralino mass and annihilation chan- 
nel. The produced neutrinos are let to interact on their way to the detector and the charged 
current interactions close to the detector where the muons are produced are simulated. The 
muons can then interact and scatter on the way to the detector. The (differential) muon 
flux at the detector is then given by summing up all these muons and weighting each muon 
by the probability that such a muon would have been created and detected, 

a cc (E„)NL,(E°) 

dct ~ 4^D2 ( b - 69 > 



where acc is the charged current cross section, Eqs. (5.26) and (5.28), N is the number 
of nucleons per cm 3 in the material surrounding the detector, is the range of a muon 
produced with energy E® and D is the distance from the source (the Sun or the center of 
the Earth) to the detector. 

This way the muon fluxes in units of m -2 annihilation -1 are obtained for the set of 



masses and annihilation channels given in Section 5.2.5. When a muon flux is needed for 
another mass, an interpolation is performed and when muon fluxes from other than these 
'fundamental' annihilation channels are needed, e.g. ZH^ , the flux is easily calculated based 
on the 'fundamental' annihilation channels. A Higgs boson will decay in flight to any of the 
particles for which the muon fluxes are calculated, e.g. bb. The Higgs bosons are let to decay 
in flight and the fluxes are obtained by integrating over the production angle of the decay 
products with respect to the Higgs momentum. The total flux from the Higgs boson is then 
obtained by summing over all the Higgs decay channels. 

In the next section, it is described how the muon fluxes for specific MSSM parameters 
are obtained. 



5.3 Muon fluxes - predictions 

We are now ready to apply the results obtained earlier in this chapter and calculate the 
expected muon fluxes at neutrino telescopes |6l| l . For other references on predicted rates in 
neutrino telescopes, see e.g. Ref. |p2| . 

The differential muon flux at a neutrino telescope is given by 

-T A A cS ^B l —— (5.40) 
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Figure 5.4: The predicted muon flux versus the neutralino mass coming from neutralino 
annihilation in a) the Earth and b) the Sun. The horizontal line is the Baksan limit |3^] . 
Only models with 0.025 < il x h 2 < 1 are shown and the muon energy threshold has been 
assumed to be 1 GeV. 




Figure 5.5: The predicted muon flux versus the relic density, £l x h 2 , for neutralino annihi- 
lation in a) the Earth and b) the Sun. The horizontal line is the Baksan limit |3£| and the 
vertical dashed lines indicate the cosmologically interesting region, 0.025 < fl x h 2 < I. The 
muon energy threshold has been assumed to be 1 GeV. 



where Ta is the annihilation rate, Eq. ( |5.2| ), A e s is the effective area of the detector in the 
direction of the source, Bi is the branching ratio for annihilation into annihilation channel 
i and d 2 <fi^/d9dE fl is the differential muon flux from annihilation channel i as obtained in 
the previous section. 

In this chapter, we have up till now only assumed that the dark matter candidate we 
consider is a WIMP, but when we want to make detailed predictions of the muon fluxes 
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Figure 5.6: The predicted muon flux from neutralino annihilation in the Earth versus the flux 
from the Sun. The horizontal and vertical lines are the Baksan limits |J8) . The dashed line, 
indicating equal rates, is shown just for convenience. Only models with 0.025 < Q x h 2 < 1 
are shown and the muon energy threshold has been assumed to be 1 GeV. 



at neutrino telescopes we have to specify what WIMP candidate we have, in our case the 
neutralino. We have performed several different scans of the supersymmetric parameter 
space as given in Table 2.2 in Section 2.£. For each given model we have checked against 
all experimental bounds given in Section 3.1 and only kept models not excluded by any 
experiment. 

As soon as the MSSM parameters are chosen, the mass of the neutralino, its composition, 
the annihilation rate and annihilation branching ratios can be calculated. The annihilation 
branching ratios Bi in Eq. ( 5.40| ) are evaluated with the same methods as those in Section 4.4 
but for the relative velocity v — since the neutralinos are highly non-relativistic when they 
annihilate in the Earth and Sun. 

We are mainly interested in models where the relic density of neutralinos is cosmologically 
interesting, i.e. 0.025 < Q x h 2 < 1, and in the figures shown here (except for Fig. we will 
only show models with a relic density in this desired range. The relic densities are calculated 
with coannihilations between neutralinos and charginos included as described in Chapter ^. 

In Fig. 5.4 we show the predicted muon fluxes coming from neutralino annihilation in 
the Earth and Sun versus neutralino mass. As seen the expected rates in neutrino telescopes 
vary over several orders of magnitude. The spread is bigger for annihilation in the Earth 
since the capture rate in the Earth only depends on the spin-independent scattering cross 
sections whereas the capture rate in the Sun also gets a contribution from the spin-dependent 
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Figure 5.7: The indirect detection rates from neutralino annihilations in a) the Earth and b) 
the Sun versus the direct detection rates in 76 Ge |pi] , pl| . The horizontal line is the Baksan 
limit |38|. Only models with 0.025 < £l x h 2 < 1 are shown and the muon energy threshold 
has been assumed to be 1 GeV. 



scattering cross sections (due mainly to hydrogen). We also show the limit on the muon 
fluxes coming from the Baksan experiment pq] and we see that present neutrino telescopes 



already have started to explore the MSSM parameter space. As will be seen in Section 5.5 
an 0(1 km 2 ) neutrino telescope can explore muon fluxes down to about 50-f00 km~ 2 yr" 1 
which at least for the Sun is a substantial fraction of the models in Fig. [TJ. Note however 
that the density of points in these figures does not have any physical meaning, they are just 
artifacts of how the scanning is performed. 

Most of our models with very high muon fluxes (i.e. those probed by Baksan) come from 
the 'light Higgs' scan in Table 2.1 where the mass of the A boson is low and hence the 
spin-independent scattering cross sections are high (which means that the capture rates are 
high). 

In Fig. |5.5| we show the expected muon fluxes versus the neutralino relic density. In this 
figure only we also show models with non-interesting relic densities. For fl x h 2 < 0.025 the 
local halo mass density p x is rescaled as p x fl x h 2 /0.025 since for these low relic densities the 
neutralinos cannot make up all of the galactic halo. The general trend of getting lower muon 
fluxes for higher £l x h 2 is clearly seen. This is due to that the relic density is approximately 
inversely proportional to the annihilation cross section (see e.g. the approximate expression 
Eq (1.6)) and due to the crossing symmetry also to the scattering cross section and hence to 
the annihilation rate. The bend over for il x h 2 < 0.025 is due to the scaling of p x explained 
above. 

In Fig. |5.6| we show the muon flux from neutralino annihilations in the Earth versus the 
flux from neutralino annihilations in the Sun. As seen, in most cases the flux from the Sun 
is higher due to that the capture rate in the Sun also gets contributions from spin-dependent 
scatterings since there is so much hydrogen in the Sun whereas the capture rate in the Earth 
only depends on the spin-independent scattering cross sections. 

In Fig. [5.7| we compare the muon fluxes in neutrino telescopes with the direct detection 
rates in 76 Ge plLplj]. Future direct detection experiments may probe event rates down to 
about 0.01 events kg -1 day -1 . We see that for a given factor of improvement in sensitivity, 
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Figure 5.8: The predicted muon flux versus the mass of the lightest Higgs boson, m H o for 
neutralino annihilation in the Sun. The horizontal line is the Baksan limit J38[ |. The upper 
part are models which can be excluded at LEP2 after 150 pb" 1 of running at 192 GeV 
and the lower part are models which cannot be excluded by LEP2. The open triangles are 
models that can be excluded due to no Higgs discovery and the open squares are models 
that can be excluded due to no chargino discovery. Only models with 0.025 < fl x h 2 < 1 are 
shown and the muon energy threshold has been assumed to be 1 GeV. 



indirect detection from the Sun generally gains more than direct detection which in turn 
usually gains more than indirect detection from the Earth. Be aware of the huge spread of 
the points though. Especially in Fig. 5.7b we see a nice complementarity between the two 
search methods. The reason for this is that spin-dependent scatterings contribute to the 
flux from the Sun but not to either the flux from the Earth or to the direct detection rates. 
There are however some direct detection materials which have spin, but the spin-dependent 
scatterings are even in these cases not contributing much to the direct rates. 

We now want to compare these rates in neutrino telescopes with the search potentials 
of LEP2. LEP2 will mainly put new constr aint s on the mass of the lightest Higgs boson, 
m H o, and on the chargino masses. In Fig. 5.8 we show the muon fluxes from neutralino 
annihilation in the Sun versus the H® mass. In the upper panel we show models that can be 
probed by LEP2 after 150 pb _1 of running at 192 GeV Jll| and in the lower panel, models 
not excludable by LEP2 are shown. For the Earth, the corresponding figure would be about 
the same but with a wider spread of the muon fluxes. Clearly, there is a very nice comple- 
mentarity between neutrino telescope search capabilities and LEP2 search capabilities. If 
we are unlucky, however, the Nature has chosen the parameters such that SUSY will escape 
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Table 5.1: The background fluxes of upward-going muons coming from atmospheric muon 
neutrinos. The units are km~ 2 yr _1 and the fluxes are within a cone of half-aperture angle 
5°, i.e. a solid angle of 0.0239 sr. 



both neutrino telescopes and LEP2. Note that LEP2 will be able to exclude all our models 
where m H o < 90 GeV. 

5.4 Backgrounds 

There are a few backgrounds for neutrino telescopes that need to be handled. The largest 
background is down-going muons produced by cosmic ray particle interactions in the Earth's 
atmosphere. These muons don't survive too far in the Earth though and can thus be reduced 
by putting the neutrino telescope deep underground. At several kilometers depth there are 
still a lot of down-going muons present, but since these are down-going and by looking for 
up-going muons (i.e. when the source is below the horizon) this background does not present 
a big problem. In the cosmic ray particle interactions in the atmosphere there are however 
also neutrinos produced [ ^2| , |63t and these will constitute the main background since they can 
not be removed by going deeper underground. By using the atmospheric neutrino fluxes in 
Ref. [ |33| and the neutrino to muon conversion formulas in Ref. pif we find the background 
muon fluxes underground given in Table 5.1. To obtain these backgrounds we have used 



the values of a and (3 valid for ice, Eq. (5.3E), and rock, Eq. (5.3£), respectively. We have 



integrated over a cone of half-apcrturc angle 5° which in many cases is a reasonable angular 
width of the signal flux. 

To discriminate the signal from the background one has to use the fact that the signal 
has an angular (and energy) dependence different from the background as will be discussed 
in the next section. From the Sun there is also a small background coming from cosmic 



ray particle interactions in the Sun's corona 65 1. This background will have a similar 
angular dependence as the signal but a different energy dependence. This background is 
of the order of 10 events/km 2 which is usually below the signals accessible with neutrino 
telescopes of order 1 km 2 being currently planned (see next section). However, in the most 
sensitive cases where the neutralino mass is high and we have a specific model to test for, 
this background should be taken into account. Since it will mostly effect the analysis when 
even bigger neutrino telescopes than 0(1 km 2 ) are built, we will for the moment neglect 
this background. 



5.5 What signal levels can be probed? 

If we now imagine having a specific neutrino telescope, how high must the signal flux be to 
be detectable? First of all it must be high enough so that we get any events and secondly it 
must be high enough to be seen above the background. Since the atmospheric background 
has both a completely different angular and energy distribution than the signal one would 
expect that it should not be too difficult to observe the signal above this background. In the 
simplest data analysis one would just look for an excess of events in an angular cone centered 
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around the Earth or Sun with a given half-aperture angle # max . The angle # max would be 
chosen large enough to include as much signal as possible and small enough to exclude as 
much background as possible. This is what is usually done with neutrino telescopes today. In 
Rcf. |36| was investigated what the optimal (9 max would be for different neutralino masses and 
compositions. Here we go one step further Q and investigate what could be done if the full 
angular and/or energy resolution is actually used to discriminate against the atmospheric 
background. Note, though, that currently designed neutrino telescopes have a very poor 
energy resolution in the energy region where the signal from neutralino annihilation is. 

In the following subsections, we describe how we calculate the expected sensitivities for 
neutrino detectors with angular and/or energy resolution and present numerical results for 
some representative models. 



5.5.1 'Single-bin' analysis 

Consider a theory which predicts an upward-muon flux <j) s (where the s stands for 'signal') 
with an angular distribution d<fi s /d0 — (psf s (6) sin#, where 9 is the angle the muon makes 
with the direction of the source of interest, and J f s (0) s'm9dd — 1 (i.e. f(9) is constant 
for an isotropic distribution). We would like to disentangle this signal from a background 
of atmospheric neutrino-induced muons which has a flux (f>f, with an angular distribution 
d(pb/d9 — <fi®fb(9) sin (9, which is nearly isotropic (at least on small angular scales). 

Consider first an experiment that can only tell that a muon has been detected with 
an angle 9 < # max and an energy > -B*' 1 , but no further information on the muon 
direction or energy is available. Then, the angular acceptance cone around the source must 
be large enough to include all (or most) of the muons produced by neutrinos from the 
source. One would therefore have some number of muons detected with an angle 9 < # max . 
For example, in their searches for energetic neutrinos from the Earth and Sun, the Baksan 
collaboration S] reports the flux of muons within an angle 6* max = 30° of the Sun or the 
center of the Earth. The Kamiokande collaboration reports the flux of muons within 
an angle varying between 6* max = 5°-30°. We will refer to this way of analyzing data as the 
'single-bin' (or '0D') approach. 

With such an experiment, the number of background events after an exposure £ (for 
example, in units of km 2 yr) is = £0° J Q e ' max fb{9) sinOdO. The number of expected 
events from the source of interest is N s = £<jfi s J max f s (9) shi9d9. A 3cr detection would 
require an excess of 3v7Vj + N s events over the number expected. Then, a 3cr excess will 
be observable only if tfi s > 3a where a = V-^fc + N s /£. 

From such a simple experiment described above where no energy or angular distributions 
are used, we can conclude that the minimal exposure required for a 3cr discovery is 

9(<fa + <fe) . . 

fmin = To (5-41) 

n 

where (f>b and <f> s are the background and signal fluxes above threshold and within the angular 



cone of acceptance # max . Note that Eq. (5.41) is only valid when the fluxes are high. This 



minimal exposure is relevant to the way, e.g., Baksan and Kamiokande have analyzed their 
data (with different values of 9 max ). 



In the following more detailed examples, Eq. (5.41) with either # max = 5° or the optimal 
$m ax (which maximizes s /<j) will be used for comparison. Note that for low masses (< 100 
GeV) the optimal # max will be higher than 5° and for high masses it will be lower. However, 
one cannot know in advance what the optimal cut will be, unless we have a specific model 
we want to test. If we don't have a specific model 5° is a reasonable choice of # max giving 
decent results both for low and high masses. In principle it is however possible to extract the 
optimal angular cut from data by varying # max and for the optimal value the signal should 
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be most clearly visible. We have not investigated further how well this method of finding 
#max would do, but it will never do better than the single-bin approach where the optimal 
angular cut is known in advance. 



5.5.2 Covariance-matrix analysis 

Now consider a slightly more sophisticated experiment which has angular and/or energy 
resolution. It is possible that we will actually be fitting for both a background and a signal 
flux of muons where the background flux is given by 



2, 



dEd6 



(E,6) = <f>° b f b {E,6), (5.42) 



which we assume to be isotropic (at least over small angular patches), fb{E, 9) = fb(E). We 
will only consider the atmospheric neutrino background resulting from cosmic-ray interac- 
tions in the Earth's atmosphere. In addition, we will want to fit data for an annihilation 
signal which generally depends on the neutralino (or any WIMP) mass m x and its compo- 
sition. We may parameterize this as 

(E, 9) = </>° [a/ hard (m x , E, 9) + (1 - a)/ soft (m x , £, 9)] , (5.43) 



dEd9 



where a parameterizes the relative contributions of a 'hard' and 'soft' annihilation spectrum. 
As a 'hard' annihilation spectrum we have used the t+t^ channel below the W mass and 
W + W~ above and as a 'soft' spectrum we have used bb. These channels represent the 
extreme hardnesses of the spectrum for any given neutralino mass. For the evaluation of 
the neutrino and muon flux for these channels we have used the method given earlier in this 
chapter. 

Therefore, we are assuming that the muon angular and/or energy distribution from both 
background and signal will be described by the set of parameters s = {4> b , (f>° s , m x , a} (one 
could also envision more parameters). We now want to ask, with what precision can we 
measure these parameters with a given experiment, assuming the true distribution is given 
by some set of parameters, s ? 

To answer this, we assume the data is binned into a number of angle/energy bins, and 
each bin i is centered on angle 9i and energy Ei with widths AEi and A0j. Therefore, for 
a given set s of parameters, the flux will be 

d ^^As) = ^(E,0;s) + ^ { E,e;s). (5.44) 



dEdO y ' ' ' dEd9 K ' ' ' dEd9" 
The probability distribution for the number of events expected in each bin is a Poisson 
distribution with mean Ni = £j^g(E i ,9 i )AE i A9 i , so it has a width a % = yfWi. 

So, suppose the true parameters are So. Then the probability distribution for observing 
an angle/energy distribution which is best fit by the parameters s is 



P(s) cx cxp 



-i(s-so) • [a] ■ (s-s ) 



(5.45) 



where the curvature matrix [a] is given approximately by 
1 dN, dNi 

&ab = t y j 



erf ds a ds b 
ds a ds b 
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Figure 5.9: The exposures needed for a 3cr discovery for different signal fluxes (indicated 
to the right in the figure in units of km~ 2 yr" 1 ) as a function of neutralino (or any WIMP) 
mass assuming perfect angular resolution but no energy resolution (and with a muon energy 
threshold of 1 GeV). The muon fluxes are evaluated for annihilation in the Earth with 
vertical background. The solid (dashed) lines correspond to hard (soft) muon spectra. The 



three signal parameters {(jr sl m Xl a\ in Eq. (|5.43|) are assumed to be unknown while the 



background flux is assumed to be known, 
yr are realistic in the near future. 



Note that only exposures less than, say, 25 km 



where the partial derivatives are evaluated at s = so, and we used cr, = \ N~i in the second 
line. In a realistic experiment, the width of the bins would be comparable to the angular 
and/or energy resolution of the experiment. In the limit of perfect angular and energy 
resolution, the sum becomes an integral, 



a ab = 48 J J dE d9 d ^ d ^ g. (5.47) 

The covariance matrix, [C] — [a] -1 gives an estimate of the standard errors that would 
be obtained from a maximum-likelihood fit to data: The standard error in measuring the 
parameter s a (after marginalizing over all the other undetermined parameters) is approxi- 
mately a a ~ C]Ja ■ If three times the standard error in the parameter <p° s is less than (fP s , for a 
given underlying model so and for a given experiment, then this model will be distinguishable 
from background at the 3cr level. 

If all of the parameters except for are fixed, then [a] is a I x 1 matrix, i.e. 1/er 2 . In 



this case, Eq. (5.47) reduces to 



To illustrate, if there were no background, Eq. (5.48) says that the statistical uncertainty 
in the number of events is the square root of the number of events, and this makes sense. 
However, if the total number of events is nonzero, then a signal has been discovered. 

In fact, if there is no background, and an event is seen, it constitutes discovery. On the 
other hand, if nothing is seen, the 95% CL upper limit to the number is 3. 
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Figure 5.10: The exposures needed for a 3cr discovery for the signal flux <jP s = 10 3 km -2 
yr _1 coming from neutralino annihilation in the Earth. The minimal exposures needed for 
a detector with neither angular nor energy resolution (0D), only angular but no energy 
resolution (ID) and both angular and energy resolution (2D ) is sh own. For the ID and 2D 
cases, results are given for all four parameters in Eqs. (5.42) and (5.43) being free (4 par), 
only the three signal flux parameters being free (3 par) and only the normalization of the 
signal flux, <fP s being free (1 par). An energy threshold of 1 GeV is used in all cases and for 
the 0D case an integration of the fluxes up to 9 max — 5° is performed. All curves are for 
hard annihilation spectra. 



5.5.3 Results 

We are now ready to perform some actual calculations using the techniques described in the 
previous subsection for the specific example of neutralino annihilation in the Sun and Earth. 
We assume that the neutrino energy spectra are either the hard or soft spectra described 
above; energy spectra from specific neutralino models should fall somewhere between these 
two extremes. Since the muon flux is proportional to the neutrino energy squared, the 
hard annihilation channels will generally be more important and hence in general the muon 
spectra will be more hard than soft. Because of the steep fall with energy of the atmospheric 
background, hard spectra generally require less exposure. In all integrations with angular 



(and energy) distribution the integration in Eq. (5.47) is performed up to 9 = 30°. For the 



atmospheric background we have used the results given in Ref. 162 



In Fig. 5.9 the minimal exposures needed to make a 3<r discovery are shown for a detector 
with perfect angular resolution but no energy resolution. The two extreme cases of soft and 
hard annihilation spectra are shown. In producing Fig. |5.9| we have assumed that only the 



signal flux parameters {(f>®,m x ,a} in Eq. ( 5.43| ) are unknown. This is not unreasonable 



since the background can be expected to be measured well by an off-source measurement. 
For annihilation in the Sun the curves are similar but less steep. We find that a neutrino 
telescope with exposures of about 1-25 km 2 yr would be able to detect signal fluxes down 
to about 50-100 km~ 2 yr" 1 . 

A comparison with what one could detect if one also has energy resolution is given in 
Fig. 5.1C where a comparison between how many parameters in Eqs. Q5.42 ) and ( |5.43 ) that 
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Figure 5.11: The muon fluxes coming from neutralino annihilation in a) the Earth and b) the 
Sun that can be discovered (at the 3<r level) with a neutrino telescope of present Amanda 
size. The exposures assumed are a) £ = 0.05 km 2 yr and b) £ = 0.03 km 2 yr. The neutralino 
annihilation spectrum is assumed to be hard. The angular cut for the 0D case has been set 
to the optimal one and only the signal parameter <j>® in Eq. (5.43) is assumed to be unknown 
in the ID case. The background flux is assumed to be a) vertical and b) horizontal as is the 
case for Amanda. The muon energy threshold has been assumed to be Ej^ 1 = 25 GeV. 



are known is also shown. We find that by having energy resolution we can gain as much as 
another factor of two in what signal fluxes we can detect. 

In case of a detector with angular resolution but no energy resolution we have also 
investigated what could be gained by varying the energy threshold and the gain is very 
small and only for neutralino masses above 100 GeV. Increasing the threshold above 10 
GeV gives no further improvements. On the other hand, large detectors (like Amanda) 
which have a threshold of tens of GeV will not lose much sensitivity either, for neutralino 
masses above 100 GeV. 



5.5.4 Example of a neutrino telescope 

Let us now consider a more specific example of a neutrino telescope. We will consider a 
neutrino telescope with a size of Amanda (at present ~ 10000 m 2 in the direction of the 
Earth and ~ 6000 m 2 in the direction of the Sun) run for 5 years, i.e. with an exposure of 



£ = 0.05 km 2 yr and £ = 0.03 km 2 yr for the Earth and Sun respectively. In Fig. 5.11 the 
muon fluxes that can be probed by such a detector is shown. We have here assumed that 
the muon energy threshold is 25 GeV and that we have an a priori model we want to test, 
i.e. we have used the optimal angular cut for the 'single-bin' analysis (0D) and we have only 



assumed that the signal flux normalization (jr s in Eq. (5.43) is unknown for the TD' case. 
We show the muon fluxes that can be probed when the neutralino annihilation spectrum is 
hard. We see that at higher masses it is for Amanda possible to probe fluxes down to a few 
hundred km~ 2 yr -1 at high masses. 



In Fig. 5.12 we show the muon fluxes that can be probed with a neutrino telescope with 
an area of 1 km 2 run for 10 years, i.e. with an exposure of £ = 10 km 2 yr. As above, 
the optimal angular cut has been used and only (jP s is assumed to be unknown. The muon 
energy threshold has again been assumed to be = 25 GeV. We see that such a detector 



can probe muon fluxes down to about 10 km~ 2 yr -1 at high masses. In Fig. 5.12b we 
see that one would gain about a factor of 1.5 by having vertical atmospheric background 
instead of horizontal. When looking towards the Sun, a neutrino telescope at the poles will 
have almost horizontal background and a neutrino telescope close to the equator will have 
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Figure 5.12: The muon fluxes coming from neutralino annihilation in a) the Earth and b) 
the Sun that can be discovered (at the 3<r level) with a neutrino telescope with an exposure 
of £ = 10 km 2 yr. The neutralino annihilation spectrum is assumed to be hard. The angular 
cut for the 0D case has been set to the optimal one and only the signal parameter <j>® in 



Eq. (5.43) is assumed to be unknown in the ID case. In a) the atmospheric background is 
vertical and in b) curves for both vertical and horizontal backgrounds are shown. The muon 
energy threshold has been assumed to be E^ 1 — 25 GeV. 



horizontal background in the beginning/end of the night and nearly vertical in the middle 
of the night (how vertical depends on the latitude and season). 

In case all three signal flux parameters in Eq. (5.43) are unknown, the curves in Figs. |5 . 1 1| 
|5.12 would be about a factor of 1.5-2 higher. If the energy threshold is higher than 25 GeV, 
the curves will be higher at low masses. 



Note that for especially the 0(1 km 2 ) telescope, as shown in Fig. 5.12, we are in the 
background dominated regime which means that if we increase the exposure by a given 
factor, the muon fluxes we can probe decrease by the square root of that factor. 



5.5.5 Discussion 

We can conclude that if want to test a specific model we only gain about 10-25% by using 
the full angular resolution compared to the 'single-bin' approach. If we don't have a specific 
model to test for, which is usually the case, we have to choose an angle # m ax in the 'single-bin' 
approach compared to which we can gain up to a factor of 2 by using the angular resolution 
as proposed here. By using the parameterization of the signal flux, Eq. (5.43) we can also 
gain some information on the neutralino mass and the hardness of the spectrum though. 

By varying the energy threshold not much more is gained, but by having energy resolution 
about a factor of 1.5-2 can be gained, slightly more at low signal fluxes and slightly less at 
high signal fluxes. 

We also note that for neutrino telescopes with a size of about 1 km 2 and only angular 
reslution, the signal fluxes we can expect to probe is in the region of 50-100 km -2 yr -1 



when all three signal parameters in Eq. (5.43) are unkonwn and about a factor of two lower 
if only the normalization is unkown. If we in addition have energy resolution we gain about a 
factor of 1.5-2 more. The signal fluxes within reach are almost an order of magnitude larger 
than the expected background coming from cosmic ray interactions in the Sun (about 10 
events km -2 yr -1 [p5[). Therefore it is quite safe to neglect this background at the present 
stage. When detectors are getting even bigger it will however be a severe limitation when 
looking for the neutrino flux from the Sun since this background is also highly directional 
as the signal. The energy dependence is quite different though, so having energy resolution 
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Figure 5.13: The projected angular distributions of neutrino-induced muons from neutralino 
(or any WIMP) annihilations in a) the Earth and b) the Sun for neutralino masses of 50 GeV 
(solid), 100 GeV (dashed), 200 GeV (dotted), 350 GeV (dash-dotted) and 750 GeV (wide 
dotted). The distributions shown are for hard annihilation channels, W + W~ for 100-750 
GeV and t + t~ for 50 GeV, and with a detector muon threshold of £^ h = 10 GeV and a 
detector (projected) angular resolution of 1.4°. 



may be beneficial in this case. 



5.6 What does a detected signal tell? 



We have in Eqs. (5.4)-(5.5) and Fig. |5.2j see n th at the width of the neutrino distribution 



x -ad ~ —^== , m x > 10 GeV. (5.49) 



depends on the neutralino mass. From Eqs. (5.4)-( p.5j ) we can derive that the root mean 
square value of the projected angular distribution given in Fig. 5.2 is 

^ 1 ^ 2£ 

V2R® y/m x /GcV 

If we could measure this width we would, since the width only depends on the neutralino 
(or any WIMP) mass, be able to determine the neutralino mass in a model-independent 
fashion. It is however not possible to measure this width directly since we measure the 
muons and not the neutrinos in neutrino telescopes. When converting the neutrino flux to 
a muon flux, we will introduce a model dependence, since we need to know the neutrino 
energy distribution which depends on the branching ratios to different annihilation channels. 
The angle the muon makes with respect to the neutrino in the charged current scattering 
decreases as the square root of the neutrino energy, i.e. we will get a smearing of the neutrino 
distribution which is smaller the harder the annihilation spectrum is and the heavier the 
neutralino is. The muon will also undergo multiple Coulomb scattering on it way to the 
detector which further smears the angular distribution. Note that both the intrinsic neutrino 
angular distribution and the charged current scattering angle tends to widen the angular 
distributions the lower the neutralino mass is. 

As we did in the previous section, we will consider the extreme cases of 'hard' and 'soft' 
annihilation spectra where we for the hard spectra use W + W~ when m x > mw and r + r _ 
when m x < mw and for the soft spectra we use the bb annihilation channel. These cases 
should represent the extreme cases for any neutralino (or any WIMP). 
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Figure 5.14: The neutralino mass versus full width half maximum #fwhm of the neutrino- 
induced muon distribution for soft and hard channels. The solid line corresponds to _E* h 
= 10 GeV and the dashed line corresponds to Ej^ = 2 GeV. The curves given are for 
neutralinos annihilating in the Earth. 



We have calculated the resulting muon angular distributions for hard and soft spectra and 
for different masses using the methods described earlier in this chapter. In Fig. 5.13 we show 
the resulting muon distributions for hard annihilation spectra, a muon energy threshold of 
i?* h = 10 GeV and a detector (projected) angular resolution of 1.4°. If we compare Fig. 5.13a 
with Fig. 5J2 we clearly see the widening of the distributions that the charged current and 
the multiple Coulomb scatterings give rise to. For neutralino annihilation in the Sun, the 
neutrino distributions are very narrow peaks at 6 V — 0° and hence the distributions in 
Fig. 5.13 reflects the smearing due to the charged current and multiple Coulomb scatterings 
only. 

Using these kind of muon distributions we can now investigate if it is possible to extract 



the neutralino mass from the width of the muon angular distributions. In Fig. 5.14 we 
show the neutralino mass versus the full width half maximum of the angular distribution 
for neutralino annihilation in the Earth. These curves are evaluated for a neutrino telescope 
with perfect angular resolution. When the resolution is more typical, 1-2°, the curves bend 
upwards at the lower-angle end thus reducing the upper limit on the mass for which the 
mass can be inferred. We find that for neutralino masses < 400 GeV, the Earth angular- 
distribution can be used to infer the neutralino mass. If the threshold is low, < 5 GeV, the 
Sun angular distribution can also be used. 

If one would observe a signal from both the Earth and Sun, one could imagine sub- 
tracting the Sun angular distribution from the Earth one and thus achieve the desired 
model-independent neutrino angular distribution. This is however not possible to do in a 
completely model-independent way since the neutrino interactions in the Sun depend on 
the neutrino spectrum and hence on model parameters. It is probably possible, however, to 
reduce the model dependence with this technique. 
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Chapter 6 

Conclusions 



We have performed a detailed evaluation of the relic density of neutralinos (in the framework 
of the MSSM) including all two-body final states at tree level and coannihilations between 
all neutralinos and charginos lighter than 2.1m x . We have found that the neutralino density 
is cosmologically interesting (0.025 < O x /i 2 < 1) for a wide range of neutralino masses and 
compositions. 

The coannihilation processes we have included are important not only for light higgsinos, 
but whenever \fi/Mi\ < 2, which includes all higgsino-like neutralinos as well as some mixed 
and gaugino-like neutralinos. In these cases coannihilations can reduce the relic density by 
typically a factor of 2-5, but sometimes even with up to a factor of 100. We have also found 
that coannihilations can increase the relic density up to a factor of 3. 

There has been claims in the literature p5, 27 that light higgsino-like neutralinos are 
never cosmologically interesting, but we found, in agreement with Ref. ]l2| ], that neutralinos 
with masses m x ~ 75 GeV and tan/3 < 2 can be cosmologically interesting. 

We have also shown that higgsino-like neutralinos with masses m x > 450 GeV can have 
Q, x h 2 > 0.025 and that the upper limit on the neutralino mass for the neutralinos not to 
overdose the Universe increases from 3 to about 7 TeV when coannihilations arc included. 
It may be noted, though, that O(TeV) massive neutralinos may appear unnatural in the 
sense that they require fine-tuning of the parameters. 

To evaluate the relic density correctly will be even more important in the near future when 



the cosmological parameters can be expected to be measured quite accurately |68|, 69 1 and 
we might want to draw some conclusion on the MSSM parameters from such measurements. 

We have also evaluated the indirect detection rates in neutrino telescopes coming from 
neutralino annihilations in the Sun and Earth. The detection rates in neutrino telescopes 
are found to, for many models, be explorable by the next generation neutrino telescopes, 
with sizes of 0(1 km 2 ) especially for neutralino annihilation in the Sun. We have also shown 
that there is a nice complementarity between which models neutrino telescopes can probe 
and which can be probed at e.g. LEP2. There are however models giving rise to very small 
detection rates in neutrino telescopes and that are not explorable by LEP2. 

We have also investigated what muon fluxes a given neutrino telescope can probe and 
we have found that neutrino telescopes with an exposure of about 1-25 km 2 yr can probe 
signal fluxes down to about 50-100 km~ 2 yr -1 if the full angular distribution of the signal is 
used. If supersymmetry is found in some other experiment and we know the relevant MSSM 
parameters, neutrino telescopes of this size will be able to search for such a specific signal 
down to about 10-50 km -2 yr . 

We have also shown that if a signal is seen, the width of the angular distribution can 
be used to infer the neutralino mass. If a signal is seen from the Earth, neutralino masses 
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m x 400 GcV can be determined. For a signal from the Sun, the mass can be inferred only 
if the muon energy threshold is small (< 5 GeV). 

We can thus at this stage conclude that very exciting times will soon come when the new 
bigger neutrino telescopes start operating. 



Appendix A 

Feynman Rules for the MSSM 



A.l Introduction 

This is a collection of Feynman rules based on the rules given in Ref. |q, j^jlq]. They are 
slightly rewritten in a format more suitable for numerical implementation and we have used 
the conventions given in Chapter ex cept for the convention on the matrix that diagonalize 
the neutralino mass matrix, Eq. ( 2.1 d| ) . 



In our actual calculations we have used the complex matrix Nij with positive mass 
eigenvalues as explained in Chapter ||. Another convention is to instead have a real matrix 
Zij in which case the mass eigenvalues can be either positive or negative. In this convention 



the neutralinos are instead of Eq. (2.11) given by 

X? = Z a B° + Z l2 W 3 + Z i3 H° + Z l4 Hl (A.l) 
The relation between the ^-matrix and the A-matrix is [|| 

Nij = V^iZij (A.2) 

where e\ is the sign of the i:th mass eigenvalue (i.e. the mass eigenvalues are e^m^ with 
mi > 0). 

The Feynman rules below are given in a form where it is very easy to choose one conven- 
tion or the other. If the convention with a real matrix Z^ is desired, one only has to drop 
the complex conjugation of Zij whenever that appears and remember that denotes the 
sign of the i:th mass eigenvalue with m 2 ; being the absolute value of the mass eigenvalue. 
If one instead wants to use the convention with a complex matrix Nij and positive mass 
eigenvalues (which is the one we have used), one just has to replace Z^ by Nj and put all 
i 1. 

Compared to Ref. ||,|(],|l6| all vertices are divided by i and the propagators are multiplied 
by i. 

In these rules, the elementary charge e > and hence the electron charge is — e. Quark 
charges are given as fractions of e and are denoted by e q , i.e. e u — +2/3 and ed = —1/3. 
Where relevant, the vertices are given as gfj^PL + g^k^R where 

Pl = ^(1-75) (A.3) 
Pr = ^(1+75) (A.4) 

are the left- and righthanded projection operators. The indices on the g's indicate the 
particles involved, for which we have used the following shorthand notation, 
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X any neutralino or chargino. 

q any quark (we only give the rules for one generation, but following Chapter ^ it is 
straightforward to extend this to three generations) , 

q any squark, 

H any Higgs boson, 

V any gauge boson. 

Note that the Feynman rules for leptons and sleptons are not given explicitly. They are 
however easily obtained from the quark and squark Feynman rules by the following substi- 
tutions 



(A.5) 



Note that vr docs not exist. 

We want to be able to evaluate different diagrams with general vertex couplings and only 
in the numerical code insert the values of the vertices for the actual particles involved. It is 
then convenient to have a convention for the order of the indices on the vertices. We have 
chosen to work with the following convention |fn| : 

• The first fermion has a bar and is hence outgoing. 



u 


-> n 


d 


-> l- 






d-L/R ~ 


Il/r 



If two Higgses are present, the order of the indices is the order in which they appear in 
the Lagrangian with the first one complex conjugated. This means for example that 
for the VHiHj vertex, gvHiHj is defined as 



£ = gvH iHj V fi Hlid^H j (A.6) 

which gives rise to the following Feynman rules 

. Hi 



Pi 



v VWW\ 9v Hi *M+»y 



pj 



or 



Pi 



, B j 



(A.7) 



v WVW\ 



-gVHiHjiPi +PjY 



Pi '* 



Hi 



(A, 



For the charged Higgses, the charge depends on the direction of the momentum of the 
Higgs. If Hi or Hj is a charged Higgs it is an H + in the first rule above and an H~ in 
the second one. For the rules given in the following sections the first Higgs is outgoing 
and the second one is ingoing and the charges are assigned accordingly. 
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G3 



• If one Higgs boson is present it is complex conjugated, i.e. if it is charged it is an H + 
moving out or an H~ moving in. 

• If three Higgs bosons are present of which two are charged, the first charged one is 
complex conjugated. 

In the following sections, the Feynman rules for all three-particle vertices are given and 
in Section A. 11 the propagators are given. 



A. 2 H-x~x vertices 
A.2.1 H^Xi 

First let us define 



gQ'L 



■i j 



2 [ z ia(gZj2 - g'Zji) 



9S"j = - [Zu(gZ j2 - g'Zjx] 



Zj3(gZ i2 - g'Zn) 
Zji{gZ i2 - g'Z a ) 



(A.9) 
(A.10) 



The vertex is §, Fig. 21] 




where 



ghnj 

n R 

gHnj 



gH 1 ij £ 3 

gHiijSi 
g 



"j cos a + S"j sin a 



ghijPL + g^ijPR 



(All) 

(A.12) 
(A.13) 
(A.14) 



ttO~0~0 

n iXi Xj 
The vertex is 



where 



9H 2 ij 
9H 2 ij 



9H 2 ij 



Fig. 21] 




gH 2 ij £ 3 



(q'-j sin a + cos a 



gku p L + gH 2ll p R 



(A.15) 

(A.16) 
(A.17) 

(A.18) 



G4 
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mm 



The vertex is |, Fig. 21] 

ff 3 ° <^ 9H 3 ij P L + 9H 3 ij P R 

X? (A. 19) 

where 

= ; '/;/. ( A - 2 °) 

5H 3 ij = -WH 3 ijSi (A.21) 

.9^ - g(0^sin/3-^cos/3) (A.22) 

A.2.2 /f x+X+ 
First let us define 

Qa = \H u ^ v n (A.23) 

Su = \J\UixV j2 (A.24) 
The vertex is @, Fig. 19] 

j?j <^ gkijPL + qh^Pr 

*t (A.25) 

where 

= -.9 [Qtj cos a + Sy sin aj (A.26) 

fl-Hjy = -.9 (Oji cos a + Sji sin a) (A. 27) 

The vertex is §, Fig. 19] 

xt 



(A.28) 



A. 2. H-X~X VERTICES 



G5 



where 



9h 2 ij 

n R 

9H 2 ij 



(qIj sin a 



S 1 *, cos a 



(Qji sin a — Sji cos a 



Hlxtxt 



The vertex is i, Fig. 19] 




& aij PL + g% aij PR 



where 



9H 3 ij 



(q*^ sin j9 + i§yCos/3 
5 ^Qji sin f3 + Sji cos /3J 



A.2.3 H'xtx] 
The vertex is @, Fig. 20] 




9h- l3 p L +9h-i 3 p R 



where 



0ff-y = -9 cos 13 



n R 

9H-ij 



-gsin/3 



Z l3 Uji - \ -(Z i2 + Zn ta,n9 w )U j2 



(A.29) 
(A.30) 



(A.31) 

(A.32) 
(A.33) 



(A.34) 

(A.35) 
(A.36) 



Given the vertex above one can derive ||70|l the following Feynman rules with changed direc- 
tion of the arrows 




0a 



" N X 



(A.37) 



GG 
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(A.38) 



(A.39) 



where C is the charge conjugation matrix and a and are spinor indices explicitly given for 
the cases of clashing arrows. When working with clashing arrows one should either use the 
spinor indices explicitly or note that if one writes down the Feynman amplitude according 
to the direction of the Dirac fermion (the chargino), the vertices with clashing arrows will 
appear as given above and if the direction is instead chosen as the one given by the neutralino 
the vertices will appear transposed. 



A. 3 V-x~X vertices 



A.3.1 W+x°xj 

First let us define 

1 

V2 



O l = 



6 R 



ZiiV* 2 + Zi2V* 1 



■^ Z i3 U j2 + Z* 2 U.ji 



(A.40) 
(A.41) 



The vertex is §, Fig. 75] 




7 M {9wii p L + &ii p n) 



(A.42) 



where 



9wij 

9wij 



90* 



(A.43) 
(A.44) 



Given the vertex above one can |70| derive the following Feynman rules with changed direc- 



A. 3. V-X-X VERTICES 



G7 



tion of the arrows 




(A.45) 




7^ (g^PL + g^PB 



(A.46) 




[7" b&ijPL+gfoiPR) C] 



(A.47) 



where the note at the end of Section A.2.S is applicable here as well 



A.3.2 Z°xtx1 



3 



First let us define 



' V * V /i - 2 Vi2V V + ^ sin2 ° w 
-U* 1 U jl -~U* 2 U j2 + 6 ij sin 2 Ow 



(A.48) 
(A.49) 



The vertex is §, Fig. 75] 




r {a^PL + g^PR) 



(A.50) 



where 



9za 

n R 

gzij 



cos 8 W 3 
cos 9w 



(A.51) 
(A.52) 
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A.3.3 Z°%% 
The vertex is §, Fig. 75] 




X? 



7^ (^Pl+^Pr) 



where 



9zij 

n R 

9zi 3 



gzij 



gzij£i£j 

~9*zij 
9 

2 cos 9w 



jA) 



A.3.4 n/xtxt 

The vertex is @, Fig. 75] 




7" fe P i+ff7« P «) 



where 



A. 4 X~q-q vertices 

Note that if we have the vertex 




9± fj PL+9f, fj P R 



then it follows 



that the vertex with reversed arrows is given by 
/ 



9f fj P L + 9^Pn 




(A.53) 

(A.54) 
(A.55) 
(A.56) 



(A.57) 

(A.58) 
(A.59) 



(A.60) 



(A.61) 



A. 4. X-Q-Q VERTICES 



G9 



A. 4.1 X + qq 

The vertex is §, Fig. 22] 



xt 



9s -PL+gf Pr 



where 



d L uj 



gm u v* 2 



V2mw sin [3 

-gUji 



Xjd R u 

The vertex is §, Fig. 22] 



(A.62) 

(A.63) 
(A.64) 



xt 



9s PL+gf .Pr 



where 



J d R u] 
'd R uj 







gm d U. 



.7 2 



\/2rriw cos (3 
The vertex is M, Fig. 22] 



(A.65) 

(A.66) 
(A.67) 



xt 



[9i Ldj PL+9l dj PR}C 



where 



Uu L dj 
n R 

i>u L dj 



gm d U. 



j- 



V2mw cos P 

-gVji 



"L 



(A.68) 

(A.69) 
(A.70) 
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XjURd 

The vertex is §, Fig. 22] 



where 



9u R dj 
n R 

yu R dj 



grriuVj 



j2 



V2mw sin /? 



(A.71) 

(A.72) 
(A.73) 



Xjdu L 

The vertex is §, Fig. 23] 



"L 



where 



9jdu L 
n R 

hfj du l 



9Vn 



gm d U. 



j 2 



\f2mw cos (3 



(A.74) 

(A.75) 
(A.76) 



The vertex is §, Fig. 23] 



Ur 



C-'lgfmnPL+gfdunPB 



where 



9jdu R 



9jdu R 



gm u V 



P- 



V2rriw sin (3 



(A.77) 

(A.78) 
(A.79) 
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A. 4. 2 x°QQ 
First let us define 



z j4 



Zji cos 9w + Zj2 sin 9y/ 
—Zji sin 9w + Zj2 cos 6y/ 



Z 



J 3 



Z 



J4 



(A.80) 
(A.81) 
(A.82) 
(A.83) 



X°jULU 

The vertex is @, Fig. 24] 




where 



7 * . 



gm u Z. 



J4 



V^mw sin /3 



~L p i p 



COS 9w 



gV2 /i . 2 \ 

^— - - e u sin 0w 2j 2 

cos 9w \ 2 / J 

[(2e„ — 1) smdwZji + cosOwZji] 



(A.84) 



(A.85) 



(A.86) 



X°u R u 

The vertex is §, Fig. 24] 



where 



qi - 



iURUJ 



M 5e u sin 6>w 
eeuZj! - 



V2 

g sin 6>w/ 



cos 6*w j2 



gm v 



y/2mw sin 



J j4 



(A.87) 

(A.88) 

(A.89) 
(A.90) 
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X°d L d 

The vertex is 



Fig. 24] 



where 



a. 



d, L dj 
R 

d L dj 



X°jd R d 

The vertex is 
d~R . 



V2m^cosp 



-V2ee d X- 1 H — - 

costly \2 

[(2e<j + 1) smO\yZji — cos 6\yZj2\ 



J cos 6% \2 ' 



cos #w 



Fig. 24] 



(A.91) 

(A.92) 

(A.93) 
(A.94) 



where 



J d R dj 



9 L d R d 3 P ^+9l d3 Pn 



V2 
V2> 



eedZ >* _ 9_ed^Jw z 



cos 9w 



J2 



e d g sin 6> w 

cos #vy 
gm d 



Z j3 



\2mw cos /3 

A. 5 H-q-q' vertices 
A. 5.1 H°qq 

These vertices can be written J(| Fig. 7 and 

q 




9HnnPL + 9§ qq PR 



(A.95) 

(A.96) 

(A.97) 
(A.98) 
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where 



9 Hi uu 


- n R 

— 9 Hmu — 


gm u sin a 
2mw sin (3 


9H 1 dd 


- n R - 

— 9n x dd — 


guid cos a 
2mw cos (3 


9 H 2 uu 


- n R 

— 9h 2 uu — 


gm u cos a 
2mw sin /3 


9H 2 dd 


- n R - 

— 9n 2 dd — 


guid sin a 
2mw cos (3 


9h 3 uu 


- -n R 

— 9h 3 uu 


igm u cot /3 
2mw 


9H 3 dd 


n R 

— 9H 3 dd ■ 


ignid tan (3 
2mw 



A.5.2 H~du 

This vertex can be written || Fig. 8] 

d 



H- 




where 



9 L h- 



n R 

9n-du 



\/2mw 
9 



V2: 



mw 



md tan (3 
m u cot (3 



H+- 




where 



L 

9n+ua 



n R 

9H+ud 



\[2 raw 
9 

\/2mw 



m u cot (3 
md tan [3 



9 H -du P L+9H-du P R 



A. 5. 3 H+ud 

From the vertex above it follows that ^ Fig. 

u 



9 H +ud P L+9H+ u d P R 



(A.100) 
(A.101) 
(A.102) 

(A.103) 

(A.104) 

(A.105) 
(A.106) 



(A.107) 

(A.108) 
(A.109) 



(A.110) 

(Alll) 
(A.112) 
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A. 6 V-q-q' vertices 

These vertices can be written || Fig. 71] 

9 



V 




where 



L 

9\Ydu 

9wdu 
9Wud 



9Wud 
L 

9Zuu 



n R 

9Zuu 



9zdd 



9zdd 



9j qq - 

9_ 

V2 



-ee 



= 



9 

V2 



2 cos 0\ 



-(1 — 2e u sin 2 



1W) 



9 ■ 2n 

e u sin Ow 



cos 9w 
9 

2 cos 9w 



(1 + 2 e<i sm 2 6 W ) 



cos 9w 



e ( i sm V\y 



A. 7 i7-scalar-scalar vertices 
A. 7.1 HHH 

These vertices can all be written |(| Fig. 9 and 10] 

H 



(A.113) 

(A.114) 
(A.115) 

(A.116) 
(A.117) 

(A.118) 
(A.119) 

(A.120) 

(A.121) 

(A.122) 



H 



m W9HHH 



where 



9HiH+h- 
9h 2 h+h- 



H 



cos(/3 — a) 



sin(/3 — a) + 



m z 



2mw cos 0w 

m z 
2mw cos 9w 



cos 2(3 cos(/3 + a) 



cos2/3sin(/3 + a) 



2>gm z 



2mw cos 0y/ 



cos 2a cos(/3 + a) 



(A.123) 

(A.124) 
(A.125) 
(A.126) 
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QH2 Hi H-2 
Q H\ H2 H% 

9HiH 3 H 3 
9H2H1H1 

9H2 H3 if 3 



3gmz 



2mw cos 9w 
gm z 

2mw cos 9w 

- cos2acos(/3 + a)] 



cos 2a sin(/3 + a) 
[2 sin 2a sin(/3 + a) 



gmz 



2mw cos 6w 
gmz 

2mw cos 8w 

+ cos 2a sin(/3 + a)] 

gm z 



cos2/?cos(/? -I- a) 
[2 sin 2a cos(/3 + a) 



2myy cos 9 



ir 



cos2/?sin(/3 + a) 



A.7.2 #gg 

These vertices can all be written J|, Fig. 11-14] 



(A.127) 

(A.128) 
(A.129) 

(A.130) 
(A.131) 



H 



mwgHqq 



where 



9h+Slu l 
9H+d R u R 

9h+S, r u l 

9 'jf + dLUR 

gH 1 u L a L 

9H 1 u r uh 

9Hiurul 
9H 1 d L d L 

9H 1 d a d a 



9 

V2 
gm u m d 



sin 2(3- 



m 2 d tan (3 + m 2 cot (3 



w 



V2m 2 



(cot (3 + tan (3) 

w 
gm d 



-{{i + A d tan/3) 
lm\ 

gm u 

V2m 2 



- (ji + A u cot (3) 



w 



gm z 



1 

mw cos 9w V 2 
■ sin a 



e u sin 2 6*w ) cos(a + (3) 



!>>>:. 



m w sin (3 ' 
gm z 
mw cos 9w 



e u sin 2 9w cos(a + /3) 



5m 2 



m w sin /3 
gm u 
2m w sin /? 



sm a 



[A u sin a — fi cos a] 



mw cos #w V. 2 



- + erf sin ) cos(a + (3) 



gm\ 



m w cos (3 
gm z 
mw cos 6*n/ 



cos a 



ed sin 2 £?iy cos(a + (3) 



(A.132) 

(A.133) 
(A.134) 
(A.135) 
(A.136) 

(A.137) 

(A.138) 
(A.139) 

(A.140) 
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9H 2 u L u L 



9H 2 u r u r 
9H 2 u r u l 

9H 2 d L d L 

9H 2 d R d R 
9H 2 d R d L 
9H 3 u L u R 

9H 3 d L d R 
9H 3 u R u L 

9H 3 d R d L 



9 m d 



m w cos (3 

gm d 
2m 2 v cos (3 

9 m z 
mw cos 6w \2 

■ cos a 



[Ad cos a — fi sin a] 

e u sin 2 6\y ) sin(a + (3) 



gml 



m w sin (3 



gm z 



mw cos 9w 
gm u 



2m 2 v sin (3 
gm z 



e u sin 2 9w sin(a + f3) 
[A u cos a + \i sin a] 



g m l 



m w sin P 



cos a 



mw cos Vw \2 
■ sin a 



+ ed sin 2 W I sin(a + (3) 



gm 2 d 



m w cos (3 



gm z 



mw cos Vw 

gm d 
2m w cos (3 



-ed sin 2 #w sin(a + (3) + 
[A^ sin a + /it cos a] 



TO 2 ^ COS /3 



sma 



Hi' 

igm 



2m 2 w 



-(A u cot/3 + /x) 



—{A d tan (3 + fi) 



2m 
igm u 
2m 2 w 



iqmd , . „ , 
(A d tan [3 + //) 



(A u cot/3 + ^) 



2to 



(A.141) 
(A.142) 

(A.143) 

(A. 144) 
(A.145) 

(A.146) 

(A.147) 

(A.148) 

(A.149) 

(A.150) 

(A.151) 

(A.152) 
(A.153) 



A. 8 V-scalar-scalar vertices 
A.8.1 VHH 

These vertices can all be written Fig. 1 and 2] 

. H out 

Pout / 



v VWWAs 



9VH out H in {Vout +PinY 



where 



gwH+Hr 

gWH!H+ 



-~sin(a-/3) 
--sin(a-/3) 



(A.154) 

(A.155) 
(A.156) 



A. 9. H-V-V VERTICES 
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9WH+H 2 — 


— — cos(ct — /3) 


( A 1 c;7"l 


9WH 2 H+ — 


--cos(a-/3) 
z 


(A.158) 


9WH+H 3 = 


i£ 
2 


(A.159) 


gWH 3 H+ = 


ig 

2 


(A.160) 


gZHiHt = 


igsm(a — (3) 
2 cos 9w 


('a 1 «i *\ 

^A.lOi J 


gZH 3 H 2 = 


igcos(a — (3) 
2 cos 9w 


(A.162) 


9ZH+H+ = 


g cos 29w 
2 cos 6*vi/ 


(A.163) 


9jH+H+ = 


— e 


(A.164) 



A.8.2 Fgg 

These vertices can all be written ||, Fig. 72] 



Pout _ 



V 



\aaaaa' 



9Vq'q(Pout + PinY 



Pin 



where 



#7 9K 9B 

gzu L u L 
gzu R u R 

9zd L d L 
9zd R d R 



-ee. 
-ee. 



'V2 
9 



(-1 + 2e u sin 2 9 W ) 



9 - -»*e 



2 cos 0^ 

t— e„sin 

cos 

— ^— (l + 2e d sin 2 6> w ) 
2 cos 6% 

e,f sin 2 6w 



cos 6*w 

A. 9 if-V-V vertices 

These vertices can all be written ||, Fig. 3] 

V 



H 




mwgHvvg 



(A.165) 

(A.166) 
(A.167) 

(A.168) 
(A.169) 
(A.170) 
(A.171) 
(A.172) 
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where 



gHiWw 


— gcos([3 — a) 


(A.174) 


9h 2 ww 


— <?sin(/3 — a) 


(A.175) 


9h 1 zz 


gm z la \ 
= cos(/3 a) 
raw cos dw 


(A.176) 


9h 2 zz 


gm z . , , 
= a sm(P a) 
raw cos #w 


(A.177) 
(A.178) 



A. 10 V-V-V vertices 

These vertices (Z°W + W~ and "/W + W~) can all be written @ 

7 

fci r^'** gvw[g a0 (ki - k 2 r + g^(k 2 - h) a 



V 

a 




q 2 — my + ie 
Internal fermion: Sp(q) - 



jq — rnj + ie 
/l + mf 



q 2 — m 2 + ie 



(A.179) 



where 

gzw+w- = gcosO w (A. 180) 

g~fW+w- = e (A.181) 

A. 11 Propagators 

First define the Mandelstam variables 

s = (P1+P2) 2 = (P3+P4) 2 (A.182) 

t = (Pi~p 3 ) 2 = (P2-P4) 2 (A. 183) 

u = (p! - PA f = ( P2 -p 3 ) 2 . (A. 184) 

We then have the following propagators j[6| for internal particles 

Internal scalar: Aniq) = - ^ (A. 185) 

q l — ra H + ie 

Internal gauge boson: D Fa[j {q) = -A v (q) ( g a/3 



(A.186) 



(A.187) 



where as usual e — Tra with T being the width of the propagating particle and ra its mass. 
The fermion propagator can in the i-channel be written 

S F (q) = [{h- fo)+m f ]A t (A.188) 



A.ll. Propagators 
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where 

A * = -T^TVe (A ' 189) 

and in the w-channel be written 

S F {q) = -[(A- in) + m f ]A u = [(fa- fa) + m f ]A u (A.190) 

where 

A u = \ r — (A.191) 

u — m,f + is 

The u-diagram can, when the incoming particles are Majorana fermions (like the neutrali- 
nos), be written in the convenient form 
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Summary of the Papers 



I. The whole chain of processes from the neutralino annihilation products in the core 
of the Sun or Earth to a muon flux at a detector is calculated with Monte Carlo 
simulation techniques. 

II. If a positive signal of WIMP annihilations in the Sun or Earth is seen at a neutrino 
telescope, it is here shown that the mass of the WIMP can be inferred from the width 
of the angular distribution. 

III. The indirect detection rates at neutrino telescopes for neutralino annihilation in the 
Sun or Earth are evaluated and also compared with direct neutralino dark matter 
searches and searches at LEP2. 

IV. The minimal exposures needed to make a discovery of WIMP annihilations in the Sun 
or Earth are calculated for different WIMP masses, compositions and signal fluxes. 

V. The relic density of neutralinos is calculated including so called coannihilation pro- 
cesses between the lightest neutralino and heavier neutralinos and charginos. 
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